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ABSTRACT 

The calculation of the molecular column density from molecular spectral (rotational or ro- 
vibrational) transition measurements is one of the most basic quantities derived from molecular spec¬ 
troscopy. Starting from first principles where we describe the basic physics behind the radiative and 
collisional excitation of molecules and the radiative transfer of their emission, we derive a general 
expression for the molecular column density. As the calculation of the molecular column density in¬ 
volves a knowledge of the molecular energy level degeneracies, rotational partition functions, dipole 
moment matrix elements, and line strengths, we include generalized derivations of these molecule- 
specific quantities. Given that approximations to the column density equation are often useful, we 
explore the optically thin, optically thick, and low-frequency limits to our derived general molecular 
column density relation. We also evaluate the limitations of the common assumption that the molecu¬ 
lar excitation temperature is constant, and address the distinction between beam- and source-averaged 
column densities. As non-LTE approaches to the calculation of molecular spectral line column density 
have become quite common, we summarize non-LTE models which calculate molecular cloud volume 
densities, kinetic temperatures, and molecular column densities. We conclude our discussion of the 
molecular column density with worked examples for C 18 0, C 1 '!), N 2 H + , NH 3 , and H 2 CO. Ancillary 
information on some subtleties involving line profile functions, conversion between integrated flux 
and brightness temperature, the calculation of the uncertainty associated with an integrated intensity 
measurement, the calculation of spectral line optical depth using hyperfine or isotopologue measure¬ 
ments, the calculation of the kinetic temperature from a symmetric molecule excitation temperature 
measurement, and relative hyperfine intensity calculations for NH 3 are presented in appendices. The 
intent of this document is to provide a reference for researchers studying astrophysical molecular 
spectroscopic measurements. 

Subject headings: ISM: molecules 


1 . INTRODUCTION 

Molecular spectral line column densities are fundamen¬ 
tal to the study of the physical and chemical conditions 
within dense molecular clouds. Chemical composition, 
molecular isomer and isotopomer ratios, molecular hy¬ 
drogen column densities, and molecular volume densi¬ 
ties are all based on a measure of the molecular column 
density. Past review s and tutorials on t h is subject, in - 
clud inglG enzell (fl 99l]h lEvans et all (U991H , lEvans •: 199 9i, 
and iGoldsmith fc Lamreil (1199911 . along with s pectro ¬ 
scopic textbooks suc h as Townes fe Schawlowl (I1975T) . 
iGordv fe Cookl (I1984T) . and Drainel ( 201111 haveprovided 
an excellent basis upon which one can derive the molec¬ 
ular spectral line column density for a variety of dense 
molecular cloud environments. The intent of this tutorial 
is to provide students and scientists a general single ref¬ 
erence point for the calculation of the molecular spectral 
line column density. 

We have organized the document such that some basic 
background information is first provided to allow a con- 
textural foundation to be lain for further calculations. 
This foundation includes a basic understanding of radia- 
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tive transfer and statistical equilibrium (Sections [ 2 ] and 
[3j), molecular degeneracy (Sections [5] and [G]), partition 
functions (Section [TJ), line strength (Sections 0 and [9]), 
and hyperfine structure (Section 1101) . The derivation of 
the most general form of the molecular spectral line col¬ 
umn density is presented in Section 0] We also explore 
a number of commonly-used practical approximations to 
the calculation of the molecular spectral line column den¬ 
sity in Section [FT] As the calculation of the molecular 
spectral line column density normally assumes constant 
excitation temperature, we address the consequences of 
this assumption when non-LTE conditions are present 
in Section m We also summarize non-LTE approaches 
to the calculation of molecular cloud volume density, ki¬ 
netic temperature, and molecular column density in Sec¬ 
tion [13] 

Molecular cloud structure and its relationship to the 
spatial resolution of a measurement system affects the 
derivation of molecular spectral line column density. To 
quantify these effects we address the distinction between 
beam- and source-averaged column density in Sectionfldl 

As worked examples are often illustrative for many as¬ 
trophysical calculations, we provide step-by-step deriva¬ 
tions of the molecular column density for several rep¬ 
resentative molecules; C 18 0, C 17 0, N 2 H + , NH 3 , and 
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H 2 CO in Section 1151 Appendices provide ancillary in¬ 
formation on some subtleties often encountered when 
calculating molecular spectral line column density, in¬ 
cluding line profile functions (Appendix [A}, conversion 
between integrated flux and brightness temperature (Ap¬ 
pendix [5]), the calculation of the uncertainty associated 
with an integrated intensity measurement (Appendix [C]), 
the calculation of spectral line optical depth using hyper- 
fine or isotopologue measurements (Appendix [D]), the 
calculation of the kinetic temperature from a symmet¬ 
ric molecule excitation temperature measurement (Ap¬ 
pendix [Ej, and relative hyperfine intensity calculations 
for NH 3 (Appendix |F|). 

2 . RADIATIVE AND COLLISIONAL EXCITATION 
OF MOLECULES 

When the energy levels of a molecule are in statistical 
equilibrium, the rate of transitions populating a given 
energy level is balanced by the rate of transitions which 
depopulate that energy level. For a molecule with mul¬ 
tiple energy levels statistical equilibrium can be written 
as: 

n i y ' Rij = y i n jRji , (1) 

3 3 


As hydrogen (H) is the most abundant atom in the in¬ 
terstellar medium, and in dense molecular clouds most 
of the H is in molecular form (H 2 ), n co iuder is usually 
assumed to be n(i? 2 )- 

At this point we can derive the Einstein relations A u i, 
B u i , and Bi u by considering only radiative excitation 
(Ciu = Cm = 0) and complete redistribution over the 
line profile (< j> u i(v ) = 4>iu(v))- Physically, this means 
that emitted and absorbed photons are completely inde¬ 
pendent. Equation |4] then becomes: 


/•OO /‘OO 

niBi u / Jv^iuiy) =n u B u i / J„4>ui(A + n u A ul 
Jo Jo 

50 /* OO 

[niBi u J v (j)iu{v)}dv= / [n u B u iJ v (j) lu (v) + n u A ul \ dv 
Jo 


’ft'iBiu'Ji'&iu —n. a B v i Ji/0i u T n u A u i. 


(5) 


For a system in thermal equilibrium, the relative level 
populations follow the Boltzmann distribution: 


^U 

ni 


9u 

= — exp 
9i 



( 6 ) 


where n, and rij are the populations of the energy levels i 
and j and Rij and Rji are transition rates between levels 
i and j. The transition rates contain contributions from: 


and the radiation field J v is described by the Planck 
Function B„(T): 


• Spontaneous radiative excitation (A,;,) 

• Stimulated radiative excitation and de-excitation 
(Rij = n, Bjj J(v)dv) 


• Collisional 

(Jl collider TJij ) 


excitation 


and 


de-excitation 


BAT) 


2 hu 3 



(7) 


Substituting Equations [5] and [7] into Equation [5] (which 
eliminates the line profile function (f>i u ) yields, after some 
rearrangement: 


where 0„(i/) is the line profile function and J v is defined 
as the integral of the specific intensity I v over the source 
of emission: 

J v = ^ j USA. (2) 

Our statistical equilibrium equation then becomes: 


y ' f TlcolliderCij Bij j J V (\)ij{y)dv J r ^ ( A, 

. j ' 0 ' j<i 

= y Tlj ^ Tl C olliderCij T Bji j J u 0ji (v)diy^j 


+ rijAji. (3) 

j>i 


For a two-level system with i defined as the lower en¬ 
ergy level l and j defined as the upper energy level u. 

= 0 and the statistical equilibrium equation 
(Equation [3| becomes: 


fll ( TlcolliderClu T Bi v 


A* OO 

/ AcpiAAdi* 

Jo 

flu ^ flcolliderCul T Bui J Ji/4*ul(l / )dl / T A u i ^ . (4) 
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which implies that: 


—Biu + B u i , 
9u 

( 8 ) 


9 iBi u — g u B u i 
2 hu 3 


A u i = 


-B 


ul • 


(9) 

( 10 ) 


For dipole emission, the spontaneous emission coefficient 
A u i can be written in terms of the dipole matrix element 
\lM u \ 2 as: 


_ 647^ 

A ul = 3^ c 3 I Wm 


( 11 ) 


3. RADIATIVE TRANSFER 

The radiative transfer equation, ignoring scattering 
processes (such as tho se involv i ng el ectr ons or dust) is 
defined as follows (see iSpitzed (1978) or iDraind (1201 it ) 
for details): 


dR T . 

— = -k u R + j v 
as 


( 12 ) 
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where 

s = Path of propagation along the line of sight 


Specific Intensity 

Absorption Coefficient 


- {niBiu ti uB u i) 

47r 

(13) 

° 2 9u~ A ( 1 gm u \ A 
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(14) 

Emission Coefficient 


hv 

—A u in u , 

47r 

(15) 


Since we generally do not know what the propagation 
path is for our measured radiation it is convenient to 
change independent variables from pathlength s to “op¬ 
tical depth” t v . which is defined as: 

dr v = n^ds, (16) 

where we use the convention adopted by iDraind (120111 ^1 
that the radiation propagates in the direction of increas¬ 
ing optical depth. Switching variables from s to r„ in 
our radiative transfer equation (Equation 1121) results in 
the following radiative transfer equation: 

dli/ — drn I u di~i/. (IT) 

where we define the Source Function S u : 


where n is the density (cm' 3 ) in the upper (u) or lower ( l ) 
energy level for a transition with frequency v. The source 
function S v is then equivalent to the Planck Function at 
temperature T ex (Equation [7]) : S v = B„(T ex ). This con¬ 
dition is sometimes referred to as “local thermodynamic 
equilibrium” (LTE). Equation ITTJl becomes: 

I v = 4(0)e _Tl/ + f exp [- (tv - t')\ BjT ex )dr'. (21) 
Jo 

If we further assume that T ex is a constant, Equation 12II 
becomes: 


L = J„(0) exp(—Ty) + BjT ex ) [1 - exp(—r„)] . (22) 

Molecular spectral line measurements involve differenc¬ 
ing the measured intensity toward a reference position 
which contains only background intensity from that mea¬ 
sured towards a molecular line source position: 


AI V = I V -IJ0) 

= IJ 0) exp(—r„) + B v (T ex ) [1 - exp(—ry)] - /„(0) 
= [BjT ex ) - BjT bg )\ [1 - exp(-T„)]. (23) 


In many cases the specific intensity I„ is replaced by 
the Rayleigh-Jeans Equivalent Temperature , which is the 
equivalent temperature of a black body at temperature 
T: 


JJT) = 


hv 

k 


eX P ( kT ) 1 


(24) 


S„ = —• (18) 

By multiplying both sides of Equation [17] by the “inte¬ 
grating factor” e Tv , we can integrate the radiative trans¬ 
fer equation from a starting point where r„ = 0 and 
If = lv(0) to find that: 

e r " (d/„ + Rdr v ) = e Tu S v dr v 

e T "I„-IjO)= f S v dr' 

Jo 

I v = hi 0)e _Tl '+ f exp [— (r„ — t')\ S v dr'. (19) 
Jo 


If we further define the Radiation Temperature Tr as 
follows: 


2 kv 2 


(25) 


(26) 


We can write our radiative transfer equation in a form 
which involves the observable Source Radiation Temper¬ 
ature Tr derived from a differencing measurement: 


T R = f [JjT ex ) - JjT bg )] [1 - exp(-T„)], (27) 


Equation [T9] is a completely general solution to the 
equation of radiative transfer (again, assuming that scat¬ 
tering is neglected). It defines the intensity measured 
by the observer (/„) as the sum of the background in¬ 
tensity (/„(())) attenuated by the interstellar medium 
(exp (—t„)) plus the integrated emission {S^dr') atten¬ 
uated by the effective absorption due to the interstellar 
medium between the point of emission and the observer 
{exp [— (tv, — t')]}. 

For an infinitely large medium the radiation field would 
be defined as blackbody: I u = £?„. We further assume 
that the medium through which the radiation is traveling 
is uniform at an excitation temperature T ex , defined as: 



1 IDraind 11201 II ) points out that ISpitzerl (119781 1 uses the opposite 
convention, that radiation propagates in the direction of decreasing 
optical depth. 


where we have introduced an extra factor / which is 
the fraction of the spatial resolution of the measurement 
fille d by the source (some times called the “filling factor”). 
See llJlich Sv, Haad ( 1976 ) for a complete derivation of the 
source radiation temperature for the case of a single an¬ 
tenna position switched measurement. 

4. COLUMN DENSITY 

In order to derive physical conditions in the interstel¬ 
lar medium it is often useful to measure the number of 
molecules per unit area along the line of sight. This 
quantity, called the “column density”, is basically a first- 
step to deriving basic physical quantities such as spatial 
density, molecular abundance, and kinetic temperature. 
Using the two-energy level system defined above, we can 
express the column density as the number of molecules 
in energy level u integrated over the pathlength ds: 

N, s j 


n u ds. 


(28) 
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Since we want to use our molecular spectral line mea¬ 
surements to calculate the molecular column density, 
which will ultimately involve the radiative transfer prop¬ 
erties of the molecular spectral line measured, we can use 
the definition of the optical depth (Equation [TCI), the def¬ 
inition of the absorption coefficient k ('Equation [TCI) , the 
Boltzmann equation for statistical equilibrium (Equa¬ 
tion [G]), the definition of the spontaneous emission co¬ 
efficient A u i (Equation HU, and our definition of the col¬ 
umn density (Equation 1281) to relate r„ to the number of 
molecules in the upper energy state N u : 


c 2 9u 
87 tv 2 gi 
c 2 9u 


A ui (j) v J ds'ni (s') ^1 - 


gin u [s ) 
g u ni(s') 


87 tv 2 gi 


AuifivNi 1 — 


T,,dv = 


87 TV 2 
8n 3 iy\m u \ 2 
3hc 

8n 3 iy\m u \ 2 

3hc 




exp 

exp 


1-1 

hv 
kT 
hu 
kT 


giN u 
g u Ni 

A u i(j)i/N u 

- 1 4> V N U 
-1 JV„, 


(29) 


where in the last step we have integrated over the line 
profile such that J 4> v = 1. Rearranging and converting 
our frequency axis to velocity (^ = in Equation l29l 
we get an expression for the column density of molecules 
in the upper transition state (N u ): 


N u 


3 h [ / hv\ 


r v dv. (30) 


At this point we have our basic equation for the number 
of molecules in the upper energy state u of our two-level 
system. In order to relate this to the total molecular col¬ 
umn density as measured by the intensity of a transition 
at frequency zq we need to relate the number of molecules 
in the upper energy level u (N u ) to the total population of 
all energy levels in the molecule N to t- Assuming detailed 
balance at a constant temperature defined by the excita¬ 
tion temperature T ex , we can relate these two quantities 
as follows: 


Ntot _ Qrot f N u \ 

A^7 ~~g7 exp \kT^J ’ 


(31) 


where we have introduced the “rotational partition func¬ 
tion” Q ro t, a quantity that represents a statistical sum 
over all rotational energy levels in the molecule (see Sec¬ 
tion 0 and g u , the degeneracy of the energy level u. 
Substituting for N u in Equation [31] the total molecular 
column density becomes: 


N tn t = 


3 h Q r 


8n 3 \fM v 


9u 


exp 



x 



r v dv. (32) 


element |/q u | 2 and associated line strength S and dipole 
moment /i (such that \g,jk\ 2 = Sg 2 ; Section[5]). For ab¬ 
sorption lines where the integral of the optical depth over 
velocity may be derived from a spectrum, Equation [35] 
may be used directly to determine the total column den¬ 
sity. For emission lines, the integral over optical depth 
is typically converted to an integrated intensity f Txdv 
(see Section HU- We derive several commonly-used ap¬ 
proximations to Equation 1321 including optically thin, 
optically thick, and the Raylcigh-Jeans approximation 
in Section HU We also evaluate the limitations of the 
common assumption that the molecular excitation tem¬ 
perature is constant (Section 1121) , summarize non-LTE 
models of the molecular column density fSection ll3l) . and 
the distinction between beam- and source-averaged col¬ 
umn density ('Section [TCI) . We then close this discussion 
of molecular column density by working through several 
example calculations (Section HU). We also discuss some 
minor issues related to the assumed line profile function 
(Appendix[Aj, the relationship between integrated fluxes 
and brightness temperatures (Appendix [Bj), and the un¬ 
certainty associated with an integrated intensity mea¬ 
surement (Appendix [C]) in the appendices. In a calcula¬ 
tion which utilizes the column density calculation formal¬ 
ism presented below, Appendix [E] describes the standard 
three-level model for low-temperature NH 3 excitation of¬ 
ten used to derive the kinetic temperature from measure¬ 
ments of the ( 1 , 1 ) and ( 2 , 2 ) inversion transitions of that 
molecule. 


5. DEGENERACIES 

For rotational molecular transitions the total degener¬ 
acy for an an energy level of a transition is given by the 
product of rotational (gj and gx) and spin ( gi ) degen¬ 
eracies: 

gu = gjgxgi■ (33) 

In the following we derive the expressions for these three 
contributions to the degeneracy of a molecular energy 
level. 


5.1. Rotational Degeneracy (gj) 

The rotational degeneracy due to the projection of the 
angular momentum on the spatial axis z, gj, exists in all 
molecules and is given by: 

gj = 2J U + 1. (34) 


5.2. K Degeneracy (gx) 

The K degeneracy (gx) describes the degeneracy asso¬ 
ciated with the internal quantum number K in symmetric 
and asymmetric top molecules due to projections of the 
total angular momentum onto a molecular axis. Because 
of the opposite symmetry of the doubly-degenerate levels 
for which K ^ 0, gx is defined as follows: 

gx = 1 for K=0 and all linear and asymmetric top 
molecules, (35) 

= 2 for 0 in symmetric top molecules. (36) 


In the following we show how to calculate the level de¬ 
generacy g u (Section [5]), the rotational partition function 
Q ro t = X) ( .9 i ex p( — It) (Section 0, the dipole matrix 


K-level doubling is due to the asymmetry in the molecule 
about the molecular axes and removes the K degeneracy 
in asymmetric top molecules. 
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5.3. Nuclear Spin Degeneracy (gi) 

The nuclear spin degeneracy gj takes account of the 
statistical weights associated with identical nuclei in a 
nonlinear molecule with symmetry (which most nonlin¬ 
ear molecules have). For a molecule with no symmetry 
or hyperfine splitting , each rotational level will have a 
nuclear spin degeneracy given by: 

ff n =H(2Ji+i)=( 2/ +if ( 37 ) 

i 

gnuclear 

9i = -, (38) 

9n 

where /, represents the spin of the *th nucleus, a is the 
number of identical nuclei, and gnuclear is given in Ta¬ 
ble □] for two of the largest classes of molecules found in 
the interstellar medium medium: those with two (C 2 „ 
symmetrj 0 ) and three (C 3 „ symmetry) identical nuclei. 
Symmetry and hyperfine splitting changes g n for all prac¬ 
tical cases. Hyperfine splitting is covered elsewhere in 
this docum ent, w hil e symmetry c onsiderations are well- 
cove red bv IGordv fe Cook! (11984T) (Chapter III.4). See 
also lTurned (|199lh for a general discussion applicable to 
the high-temperature limit. Keep in mind that if you 
are only interested in studying one symmetry state (i.e. 
the para species) of a molecule, gi = 1. In the following 
we list some examples of g nuc iear calculations for several 
molecules. 


5.3.1. H 2 CO and c-C 3 H 2 

Formaldehyde (H 2 CO) is a (slightly; k = —0.96; see 
Equation 05]) prolate (see Section 0] asymmetric top 
molecule, while Cyclopropenylidene (C-C 3 H 2 ) is a cyclic 
(slightly; k = +0.69) oblate (see Section [ 6 | asym¬ 
metric top molecule. Both have two opposing identi¬ 
cal H (spin=i) nuclei. The coordinate wavefunction is 
symmetric/asymmetric for K_i even/odd, respectively. 
Therefore, from the two identical spin nuclei cases in 
Table 0] 

gnuclear — (27 + 1)/ 1 foi K_i even, (39) 

= (2I + 1)(/ + 1) = 3 for K _1 odd. (40) 

5.3.2. NH 3 and CH 3 CN 

Ammonia (NH 3 ) and Acetonitrile (CH 3 CN) are sym¬ 
metric top molecules with three opposing identical H 
(spin=i) nuclei. Therefore, from the three identical spin 
\ nuclei cases in Table 0] 

gnuclear = | (2Z + 1)(4/ 2 + 4/ + 3) = 4 for K=3n(41) 
= i(2/+1)(4/ 2 + 4J) = 2 for K^3n. (42) 

O 

5.3.3. c-C 3 HandS0 2 


2 The number of symmetry states for a molecule are deter¬ 
mined by the number of configurations within which the wave- 
function of the molecule is unchanged with a rotation of n 
about a symmetry axis and a reflection of n through that 
symmetry plane. For a nice tutorial on molecular symmetry 
see Stefan Immel’s extensive discussion at http://csi.chemie.tu- 
darmstadt. de /ak/immel/tutorials /symmetry/index7. html 


Cyclopropynylidyne (C-C 3 H) is an oblate (k = +0.17) 
asymmetric top molecule with two opposing identical 
carbon (spin= 0 ) nuclei, while S0 2 is a prolate (k = 
—0.94) asymmetric top molecule with two opposing iden¬ 
tical oxygen (spin=0) nuclei. For both molecules the co¬ 
ordinate wavefunction is symmetric/asymmetric for K_i 
even/odd, respectively. Therefore, from the two identical 
spin 0 nuclei cases in Table 0] 

gnuclear = (2 1 + 1)1 = 0 for K_i even, (43) 

= (2/ + 1)(/ + 1) = 1 for K_i odd. (44) 

This indicates that half of the levels are missing (those 
for which K_i is even). 


6 . SYMMETRY CONSIDERATIONS FOR 
ASYMMETRIC ROTOR MOLECULES 

The symmetry of the total wavefunction ip for a given 
rotational transition is determined by the product of the 
coordinate wavefunction ip e ip v tp r and the nuclear spin 
wavefunction ip n . It is common to refer to a symmet¬ 
ric nuclear spin state as “ortho” and an anti-symmetric 
nuclear spin state as “para”. These wavefunctions are of 
two types: Fermions and Bosons. Table 0] lists the sym¬ 
metries for the various wavefunctions in both cases for 
exchange of two identical nuclei. Since an asym¬ 

metric top can be thought of as belonging to one of two 
limiting cases, prolate or oblate symmetric, we need to 
consider these two cases in the context of the coordinate 
wavefunction ip e tp v xp r . 

Limiting Prolate:: We consider the symmetry of the 
coordinate wavefunctions with respect to rotation 
of 180° about the axis of least moment of inertia. 
Since the coordinate wavefunction ip e ip v 'ip r depends 
on this rotation angle £ as exp (±iA"_i£), it is sym¬ 
metric when K-i is even and antisymmetric when 
K_ i is odd. H 2 CO is a limiting prolate asymmetric 
top molecule (though it is only slightly asymmet¬ 
ric) . 


Limiting Oblate We consider the symmetry of the co¬ 
ordinate wavefunctions with respect to rotation of 
180° about the axis of greatest moment of inertia. 
Since the coordinate wavefunction ip e ip v ip r depends 
on this rotation angle £ as exp (±*/\ + i£), it is sym¬ 
metric when K + \ is even and antisymmetric when 
K + 1 is odd. NH 2 D is a limiting oblate asymmetric 
top molecule. 


The level of asymmetry in molecules is of ten describ ed 
in terms of Ray’s asymmetry parameter n (I Ravi 11933) : 


K = 


2 B-A-C 
A-C 


(45) 


where A, B, and C are constants which represent the 
three principle m oments of i nertia for an asymmetric ro¬ 
tor molecule (see IGordv fe Cookl (11984T) for more infor¬ 
mation), usually expressed in megahertz (MHz). The 
limits to k are: 


• B = A:k = +1 and the molecule is an oblate sym¬ 
metric rotor. 


• B = C:re = —1 and the molecule is a prolate sym¬ 
metric rotor. 
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Table 1 

Nuclear Statistical Weight Factors for 02v and C 3 V Molecules 3, 


Identical Nuclei (cr) b 

Spin 

J 

K c 

9nuclear 

2 

1 3 

2 ’ 2 ’ ' ’ ’ 

Any 

Even 

(2 1 + 1)1 

2 

Any 

Odd 

(2J + 1)(J + 1) 

2 

0 . 1 , 2 ,... 

Any 

Odd 

(2/ + l)(J + l) 

2 

0 , 1 , 2 ,... 

Any 

Even 

(2 1 + 1)1 

3 

Any 

Even or Odd 

3n, ^ 0 

|(2/ + l)(4/ 2 + 4/ + 3) 
|(2/ + l)(4/ 2 +4/) 

3 

Any 

Even or Odd 

7 ^ 3n 


Derived from fGordv Sz Cookl (119841) , Table 3.2 and 3.3. 


U. C I r a 1 ' 


b or = 3- _ 

— (2/+l) CT 

c Where n is an integer. 


Table 2 

Eigenfunction Symmetries for Exchange of Two Identical Nuclei 11 



Spin (I) 


Wavefunction b 



Statistics 

Total ip ) 

Coordinate (ipeipv'ipr) 

Spin (ip n ) 

Qnuclear 

Fermi 

1 3 

2 , 2 » • • • 

A 

s 

A 

(2 1 + 1)1 

Fermi 

A 

A 

S 

( 2 /+l)(/+l) 

Bose 

0 , 1 , 2 ,... 

S 

S 

S 

( 2 / + l)(J + l) 

Bose 

0 , 1 , 2 ,... 

S 

A 

A 

(21 + 1)1 


a Frorn tGordv fc Cookl 119841 '). Table 3.2. 

b Key: A = Asymmetric; S = Symmetric. For nuclear spin states asymmetric = “para” 
while symmetric = “ortho”. 


Taking a more general example, the rotational angular 
momentum constants for H 2 CO are A = 281970.37MHz, 
B = 38835.42558MHz, and C = 34005.73031MHz. 
Equation H5l then yields k ~ —0.961, which means that 
H 2 CO is nearly a prolate symmetric rotor. 

7. ROTATIONAL PARTITION FUNCTIONS ( Q ROT ) 

For a parcel of gas that exchanges energy with the am¬ 
bient medium, statistical mechanics states that the par¬ 
tition function Q which describes the relative population 
of states in the gas is given by: 

Q = Y,3 ie ' K ( 46 ) 

Following iGordv fc Cook! (11984H (chapter 3, section 3), 
the partition function for molecules in a gaseous state 
is a function of the electronic, vibrational, rotational, 
and nuclear spin states of the molecule. Assuming that 
there are no interactions between these states, the total 
partition function for the molecule can be expressed as 
the product of the partition functions of these four types 
of energy states: 

Q = QeQvQrQn- (47) 

The ground electronic state of most molecules observed 
in the interstellar medium are stable, 1 £ electronic 
states, although there are a few notable exceptions (i.e. 
OH 2 n, SO 3 £, CN 2 £, etc.). For these molecules 
Q e = 1. For simplicity we will also assume that the 
molecules are in their ground vibrational state ( Q v = 1). 
This leaves us with rotational and nuclear partition func¬ 
tions comprising the total molecular partition function, 


which we can write as: 

rot = QrQn 

= 9J9Kgiex p > ( 4 8) 

J,K,I ' ' 

where t he d egenerac ies ff j, gx, a nd g j are described in 
Section 15.11 Section El Section 15.31 respectively. See 
iTurned (1199111 for a nice general discussion listing expres¬ 
sions for Q ro t in the high-temperature limit for a variety 
of molecules. In the following we derive the rotational 
partition function Q ro t for linear, symmetric, and asym¬ 
metric rotor molecules. 

7.1. Linear Molecule Rotational Partition Function 
For linear molecules: 

• gj = 2 J + 1 (Sect ion HTT1) 

• gx = 1 (Section [OJ) 

• < 7 / = 1 (since linear molecules with allowed electric 
dipole transitions are polar and have no center of 
symmetry) 

which implies that Equation l48l becomes: 

Qrot = ( 2J + 1 ) ex P ' ( 49 ) 

The energy levels for a linear molecule can be described 
by a multi-term exp ansion as a function of J(J + 1) 
(j.Tennings et alJfi987 ): 

Ej = h(B 0 J{J + 1) - D 0 J 2 (J + l) 2 + Ho J 3 (J + l) 3 
— + 1)^ + + l) 3 H - ...), (50) 
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Figure 1. Rotational partition function calculations for CO u sing 
the lowest 51 levels of the molecule. Shown are Qrot (Eqiiat ion l-191 
vertical scale on left) and the p erce ntage differences bet ween Qrot 
and the expansions of Equation 1191 provided by Equation 1521 (using 
terms to third-order in ) and Equation 1531 (vertical scale on 
right). 


where B 0 is the rigid rotor rotation constant and D 0 , H 0 , 
L 0 , and M 0 are the first- through fourth-order centrifugal 
distortion constants for the molecule, respectively, all in 
megahertz (MHz). Using the rigid rotor approximation 
to the level energies, thus ignoring all terms other than 
those linear in J(J+ 1 ), Equation IHUl becomes: 

Ej = hB 0 J{J + 1). (51) 

This allows us to approximate Q ro t for diatomic linear 
molecules as follows: 


Qrot ~ ^2 (2 J + 1 ) exp ( - 


j -o 


hB 0 J(J + 1) 
~kT 


,, kT 11 fhBo\ fhB 0 \ 

~hB 0 + 3 + 15 V kT J + 315 V kT J 


+ 


1 

315 



(52) 


(fro m IGordv &; Cookl (|1984l ) Chapter 3, Equation 3.64 
and iMcDowelll (Il988f) , where application of the Eulcr- 
MacLaurin summation formula has been used). This ap¬ 
proximate form is good to < 1% for T > 2K (Figure [lj. 
This level of accuracy is maintained even if one includes 
only the first two terms in Equation 1521 
An alternate appro ximation for linear polyatomic 
molecules is derived by IMcDowelll (119881 1 : 


Q 


rot — 


kT 

hB 0 


exp 



(53) 


which is reported to be good to 0 . 01 % for < 0.2 
(T > 13.7 K for CO) and is good to better than 1% for 
T > 3.5 K (Figure [TJ) . Note that Equation 1531 reduces to 
Equation [52] when expanded using a Taylor Series. 


7.2. Symmetric and Slightly-Asymmetric Rotor 
Molecule Rotational Partition Function 

For symmetric and asymmetric rotor molecules: 

• gj = 2 J + 1 (Section [OJ 

• gx = 1 for K = 0 and 2 for K ^ 0 in symmetric 
rotors (Section [5T2li 

• gx = 1 for all K in asymmetric rotors 

• gi = 9 £T+ly (See Table [TJ 

which implies that Equation l48l becomes: 

Qrot = 9 Kgi(?J + 1 ) exp (--Tp) ■ ( 54 ) 

j=ox=-j ' ' ' 

Like the energy levels for a linear molecule, the en¬ 
ergy levels for a symmetric or slightly-asymmetric [3 rotor 
molecule can be described by a multi-term expansion as 
a function of J(J + 1): 

Ej K = h(B 0 J(J + 1) + s 0 K 2 + DjJ 2 (J + l) 2 

+ D jk J(J + 1 )K 2 + D k K A + Hj kk J(J + 1 )AT 4 

+ H Jjk J 2 (J + 1 ) 2 K 2 + H 3& J 3 {J + l) 3 

+ H k6 K 6 + (55) 

where so = Aq — Bq for a prolate symmetric rotor 
molecule and so = Co — Bq for an oblate symmetric 
rotor, and the other constants represent various terms 
in the centrifugal distortion of the molecule. All con¬ 
stants are in megahertz (MHz). For rigid symmetric or 
slightly-asymmetric rotor molecules, using the rigid rotor 
approximation to the level energies: 

Ejx = h(B 0 J(J + l) + s 0 K 2 ). (56) 

From IMcDowelll (|199f)D we can then approximate Q ro t 
for a symmetric or slightly-asymmetric rotor molecule as 
follows: 


„ Vmn 

Qrot — - exp 


/iHq( 4 —m)\ / kT 


12 kT 


hBn 


3/2 


1 + 


1 / hBo(l — m) 


90 


kT 


, (57) 


where 


m =—^- for a prolate symmetric rotor molecule, 
A o 

=—)- for an oblate symmetric rotor molecule, 
t-o 


B 2 0 

~A 0 Cn 


for a slightly-asymmetric rotor molecule 


(see IH erzberd (|1945h ). 


3 Slightly-asymmetric rotor molecules are defined such that Bo ~ 
Co- In Equation 1551 replace B with y/B qCq- 






































Mangum & Shirley 



100 


80 


60 g 

o 

<5 

CL 


40 IS 

O 


20 


0 


Figure 2. Rotational partition function calculations for NH 3 us¬ 
ing the lowest 51 levels of t he m olecule. Shown are the exact sum¬ 
mation for Q rot (E quat ion I54|l and the approximate forms given 
by Equations I57l and l58l (vertical scale on left). Shown also are the 
fractional percentage differences (vertical scale on right), given by 
100 * Qrot (exact) -Q r ot (approx) q £ ^ese two approximations rela- 

Qrot (exact ) ’ 

tive to Qrot (Equation l54)l . 


If we expand the exponential and take only up to first 
order terms in the expansion in Equation 1571 


Qrot — 


IWl'K 

a 

/rmr 
a 


1 + AB 0 (4- to) + 


kT 

hBa 

kT 

hBo 


12 kT 
3/2 

1/2 


( kT 


V hBn 


3/2 


(58) 


which is the eq uation for symmetric r otor partition func¬ 
tions quoted bv lGordv fc Cookl ()1984f l (Chapter 3, Equa¬ 
tions 3.68 and 3.69). Figure [2] compares Q rot calculated 
using Equation [54] and the approximate forms given by 
Equation [57] and Equation [5§] for NH 3 . In this exam¬ 
ple Equation [57] is good to < 20% for Tk > 10 K and 
< 0.25% for Tk > 50 K, while Equation [55] is much less 
accurate, good to < 40% for Tk > 10 K and < 6 % for 
T k > 50 K. 


8 . DIPOLE MOMENT MATRIX ELEMENTS 
(l/rjifl 2 ) AND LINE STRENGTHS (5) 

The following discus sion is derived from the excel¬ 
lent discussion given in iGordv fe Cook! (11984D . Chapter 
II. 6 . A detailed discussion of line streng ths for diatomic 
molecules can be found in iTatunJ (11986I L Spectral tran¬ 
sitions are induced by interaction of the electric or mag¬ 
netic components of the radiation field in space with the 
electric or magnetic dipole components fixed in the ro¬ 
tating molecule. The strength of this interaction is called 
the line strength S. The matrix elements of the dipole 
moment with reference to the space-fixed axes (X,Y,Z) 
for the rotational eigenfunctions i[) r can be written as 


follows: 


/ = £ Mg / '^r^Fgi’rdr, (59) 

J 9 J 

where 4 > Fg is the direction cosine between the space-fixed 
axes F=(X,Y,Z) and the molecule-fixed axes g=(x,y,z). 
The matrix elements required to calculate line strengths 
for linear and symmetric top molecules are known and 
can be evaluated in a straightforward manner, but these 
calculations are rather tedious because of the complex 
form of the eigenfunction. Using commutation rules be¬ 
tween the angula r momentum oper ators and the direc¬ 
tion cosines T . iGross et al.l (!l944f ) derive the nonvan¬ 
ishing direction cosine matrix elements in the symmetric 
top representation (J,K,M): 

(J, K, M\$Fg\J', K', M') = {J\$Fg\J') 

x (J,K\<S> Fg \J',K’)(J,M\$ Fg \J',M'). (60) 

where M is the projection of the total angular momen¬ 
tum J on a fixed space axis. There are 2J+1 projections 
of M allowed. The dipole moment matrix element |/q „| 2 
can then be written as: 

M 2 = E ^2\(j, k , m \tf\j', k ', m ')\ 2 , ( 6 i) 

F=X,Y,Z M’ 

where the sum over g = x, y, z is contained in the expres¬ 
sion for g F (Equation I551l. Table [3] lists the direction co¬ 
sine matrix element factors in Equation [221 f° r symmetric 
rotor and linear molecules. In the following we give ex¬ 
amples of the use of the matrix elements in line strength 
calculations 


9. LINEAR AND SYMMETRIC ROTOR LINE 
STRENGTHS 

For all linear and most symmetric top molecules, the 
permanent dipole moment of the molecule lies completely 
along the axis of symmetry of the molecule (fj, = n z ). 
This general rule is only violated for the extremely-rare 
“accidentally symmetric top” molecule (where I x = I y ). 
For all practical cases, then, Equation I6T1 becomes: 

lw„| 2 =m 2 £ £|(j,il,m|$ Fs |j',il',m , )| 2 

F—X,Y,Z M' 

= H 2 S, (62) 

where we have defined the line strength S as: 

£ £|(J,iL,M|4> Fg |J',IL',M')| 2 . (63) 

F=X,Y,Z M ' 


9.1. (J, K) —> (J — 1, K) Transitions 

Using the matrix element terms listed in the fourth 
column of Table [3] with the definition of the direction 
cosine matrix elements (Equation 1601) . we can write the 
terms which make-up Equation [62] for the case (J,K) —> 
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Table 3 

Direction Cosine Matrix Element Factors 3, for Linear* 3 and Symmetric Top Molecules 




J' Value 


Matrix Element Term 

J + 1 

J 


{J\*Fg\J') 

{J,K\<& Fz \J' ,K) 

(J,K\<f>Fy\J',K ± 1) = ^i(J,K\<f> Fx \J',K±l) 

(J,M\$ Zg \J',M) 

(J,M\<t> Yg \J’,M±l) = ±i{J,M\$ Xg \J',M±l) 

|4(J + 1) [(2J + l)(2J + 3)]i} _1 

2[(J + f)2_iy2]^ 

^[(J ± K + 1)(J ± K + 2)]? 

2 [(J + l) 2 - M 2 ] 3 

=F [(J ± M + 1)( J ± M + 2)] i 

[4J(J + 1)] _1 

2 K 

[J(J + 1) - K(K + 1)]3 

2 M 

[J(J + 1) - M(M+ 1)]3 

[4J(4J 2 -4- l)i] 

-2(J 2 - ^2)5 

T[(J + M(JT A - 1)]3 

—2(J 2 - M 2 )3 

T[(J + M)(Jq=M- 1)]3 


3 Derived from IGordv k, Cookl 11984 1. Table 2.1, which is itself derived from ICross et all 11944 1. 
b For linear molecules, set K=0 in the terms listed. 


as follows: 


(Section [SJ : 


|W «| 2 = M 2 


(J 2 - AT 2 ) 1 / 2 


J(4J 2 - 1)V2 _ 


i r 


r at 2 1 

i (J 2 - M 2 ) 1/2 

s = 

4J 2 (J+I) 2 j 


AM 1 


+ ( 1 r 1 ) [(J=FM)(J ™ _1)]1/2 }- (64) 

Applying these terms to the dipole moment matrix ele¬ 
ment (Equation [61] which simply entails squaring each 
of the three terms in Equation [5?] and expanding the ± 
terms) and using the definition of |/q u | 2 (Equation [63]) : 


' (J 2 - AT 2 ) ' 
J 2 (4J 2 - 1) 



M 2 ) 


+ \[{J~ - M - 1) + (J + M){J + M- 1)] . 

(65) 


+2 [J(J + 1) - M(M +1) + J(J + 1) - M(M - 1)] . 

( 68 ) 

Reducing Equation [05] results in the following for a sym¬ 
metric top transition ( J,K ) —A ( J,K ): 

5 = jJ+ i) for ( J ’ K ) ( J ’ K )‘ ( 69 ) 

9.3. (J,K)—t(J + l,K) Transitions 

Using the matrix element terms listed in the second 
column of Table [3] with the definition of the direction 
cosine matrix elements (Equation 1601) , we can write the 
terms which make-up Equation [551 for the case (J, AT ) —> 
(J + l,Atl as follows: 


Reducing Equation [65] results in the following for a sym¬ 
metric top transition ( J,K ) —A (J — 1,AT): 

5 = /(2J-fl) fOT ( ' 1 L K )‘ ( 66 ) 

To derive the equation for a linear molecule transition 
J —> J — 1, simply set AT = 0 in Equation [66] 

9.2. (J, AT) —> (J, AT) Transitions 


2 = 2 [i J + l) 2 ~ A ' 2 ] 1/2 

M (J+1)[(2J+1)(2J+3)] 1/2 

x | [(J + l) 2 - M 2 ] 1/2 

- (Hr) [( J±M + 1 )( J±M + 2 )] 1/2 |- (70) 


Using the matrix element terms listed in the fourth 
column of Table [3] with the definition of the direction 
cosine matrix elements fEuuation l60l) . we can write the 
terms which make-up Equation [55] for the case (J, AT) —> 
(J,kB as follows: 


I Abu I — A 

*{ 


2 AT 

4J( J + 1) 

2M ± 2 [J( J + 1) - M(M ± 1)] 


1/2 


}■ 


(67) 


Applying these terms to the dipole moment matrix el¬ 
ement ('Equation IQTT) and using the definition of |/q u | 2 


4 Also referred to as “P-branch transitions”. 

5 Also referred to as “Q-branch transitions”. 


Applying these terms to the dipole moment matrix ele¬ 
ment (Equation [61] which simply entails squaring each 
of the three terms in Equation [70] and expanding the ± 
terms) and using the definition of |/i; u | 2 (Equation [63]) : 


S = 


[(J + l) 2 ^ 2 ] 


1/2 


(J + l) [(2J + 1)(2J + 3)] 1/2 
1 


(J+l) 2 - M 2 


+ 


(J + M + l)(J+M + 2) 


+ (J — M + 1)(J — M + 2) 


>. (71) 


Also referred to as “R-branch transitions”. 
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Reducing Equation [7l] results in the following for a sym¬ 
metric top transition (J, AT) —> (J + 1,K): 

s = WTM7TTj la{J - K ^ {J + 1 - K) - (72) 

To derive the equation for a linear molecule transition 
J —> J + 1, simply set K = 0 in Equation [72] 


10 . HYPERFINE STRUCTURE AND RELATIVE 
INTENSITIES 

Hyperfine splitting in molecular spectra occurs when a 
nucleus in the molecule has non-zero nuclear spin. Hy¬ 
perfine splitting is commonly observed in molecules that 
contain 4 H (I = 1/2), 2 D (I = 1), 13 C (I = 1/2), 14 N (I 
= 1), 15 N (I = 1/2), 17 0 (I = 5/2), and 33 S (I = 3/2) 
nuclei. For a molecule with a single coupling nucleus, the 
vector addition of angular momenta is used to define a 
new hyperfine quantum number, F, such that F = I + J 
is the total molecular angular momentum, where I and 
J are the nuclear spin and rotational angular momentum 
of the molecule respectively (I and J precess around F). 
The allowed scalar values of F are given by 


F = J + I,J + I-1,...,\J-I\ . (73) 

Selection rules stipulate that A F = 0, ±1 (but F = 0 
cannot go to 0). For example, the 14 N nucleus in HCN 
is responsible for strong hyperfine coupling because the 
spin of 14 N is I = 1. Equation 1751 shows that the J = 0 
level remains unsplit (F = 1), the J = 1 level splits into a 
triplet (F = 2, 1,0), and the J = 2 level also splits into a 
triplet (F = 3, 2, 1). The selection rule allows 3 hyperfine 
transitions from J = 1 — 0 and six hyperfine transitions 
from J = 2—1. Since 1 = 1 and the Clebsch-Gordon 
series (Equation 1731) terminates with F = |J — I|, all of 
the levels with J > 2 will also be split into triplets. 

The relative strengths ca n be calcu l ated by usi ng irre¬ 
ducible tensor methods fsee lGordv fe Cookl (11984T) Chap¬ 
ter 15). We define the relative strength such that the 
sum of the relative strength, Ri, of all transitions from 
F' —>• F for a given J' —> J are equal to one: 

Y^Ri(J'F'^ JF) = 1. (74) 

F'F 

The relative line strengths are calculated in terms of a 
6 — j symbol, 


Ri{J'F' —> JF ) 


(2E + 1)(2F' + 1) \I F' J' 
(2/ + 1) \l J F 


2 

(75) 


With the aid of 6 — j Tables (lEdmondslll960lf1 . and the 
property that 6 — j symbols are invariant with permuta¬ 
tion of the columns, we find the appropriate 6 — j symbol 


7 Many online calculation tools are available that will calculate 
6 j symbols. For example, see http://www.svengato.com/sixj.html 


for each transition: 


fa b c \ 

c- 1 6 - If ~ 

s(s + l)(s — 2 a — l)(s — 2 a) 

(26 - 1)25(26 + l)(2c - l)2c(2c + 1) 

fa b c\ 

c — 1 bf ~ 

2 (s + l)(s — 2 a)(s — 26)(s — 2 c + 1 ) 

26(26 + 1)(26 + 2 )( 2 c - l) 2 c( 2 c + 1 ) 

fa b c 1 2 _ 

|l c - 1 6 + 1 J — 

(s - 26 - l)(s - 26)(s - 2 c+ l)(s - 2 c + 2 ) 

(26 + 1)(26 + 2)(26 + 3)(2c - l)2c(2c + 1) 

fa 6 cl _ 

l 1 c 6 / - 

[26(6 + 1 ) + 2 c(c + 1 ) — 2 a(a + l )] 2 
26(26 + 1)(26 + 2 ) 2 c( 2 c + l)( 2 c + 2 ) ’ 

(76) 

where s = a + b + c. For example, in the case of the 
J = 1 — OF = 2—1 transition of HC 14 N, the first 6 — j 
symbol in Equation I7S1 should be used with a = I = 1, 
6 = F' = 2, and c = J' = 1. 

The formalism above may be generalized when mul¬ 
tiple coupling nuclei are present. By extension, an ar¬ 
bitrary number, n, of coupling nuclei may be included 
in the nested vector addition of angular momenta (i.e. 

Fi = J + I \, F ‘2 = F-\ + I 2 , • • • , F = F n _i + I n ) and 
the relative intensities can be calculated by multiplying 
the relative strength equation n times with appropriate 
quantum numbers for each step in the nested vector ad¬ 
dition. For example, in the case of two coupling nuclei 
with unequal strength, such as the outer (Ji) and inner 
(/ 2 ) 14 N of N 2 H + , the nested coupling scheme F A = J+Ii 
and F = F± + / 2 is used. The relative strength of the 
hyperfine transitions are then given by 


Ri(J'F[F' JF X F) = 

(2 F[ + l)(2Fi + 1)(2F' + 1)(2F + 1) 
( 2 Ii + l)( 2/ 2 + 1 ) 



While we have specifically discussed hyperfine splitting 
from nuclei with unequal coupling in this section, the ad¬ 
dition of vector angular momenta in quantum mechan¬ 
ics is a general problem and the formalism above may 
be extended to problems with equal coupling nuclei (i.e. 
H 2 CO) or for fine structure splitting (i.e. in molecules 
with electronic states with total electronic spin S > 0) 
with an appropriate mapping of the quantum numbers. 

We end this section by noting an important caveat 
with regards the interpretation of the relative inten¬ 
sities of hyperfine transitions. Anomalies between 
observed hyperfine transition intensities and those 
predicted by the quantum mechanics described in 
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this s ec tion have been observ e d (see iKwan fe Scoville 


1973), 

Guilloteau & Baudry C1981D, Walmslev et al. 

1982)). 

Stutzki et al. d 19841). iStutzki & Winnewissei 


( 198511 . iBachiller et, al. fl l 99 7li . iDaniel et alJ (I2006H . and 
iHilv-Blant, et al.f ( 201C ) for examples). These anomalies 


are likely due to spectral line overlap between hyperfine 
transitions and can result in non-LTE hyperfine ratios. 


11.2. Optically Thin with Rayleigh-Jeans 
Approximation 

Assume that r -C 1 and hv -C kT ex . This reduces 
the term in [ ] in Equation 1571 to Jj/(T) to T, and 
reduces the radiative transfer equation to that derived in 
Equation 1751 Equation 1521 then reduces to 


11 . APPROXIMATIONS TO THE COLUMN 
DENSITY EQUATION 

In the following we derive several commonly-used ap¬ 
proximations to the column density Equation [32j Opti¬ 
cally thin (r <C 1), optically thin when in the Rayleigh- 
Jeans limit (hv -C kT ex ), and optically thin when in the 
Rayleigh-Jeans limit and negligible background emission 
(T bg < T ex ). 

11.1. Optically Thin Approximation 

If we assume that < 1, the radiative transfer equa¬ 
tion ('Equation [271) becomes: 

T R = f[J(T ex )- J(T bg )\T, (78) 

and the column density equation fEauation l571) becomes: 


7 \jthin _ 

iy tot ~ 


exp 


\ kT ex J 


f 3 h \ / Qrot \ _ 

V8t T 3 Sfi 2 Ri) \gjgK9i) exp (Jg-) 

V Klex I 

f T R dv 


- 1 


/ (MT ex ) - j v {T hg )y 


(79) 


The J u (T ex ) — J u (T bg ) term may be removed from the 
integral because the Planck function does not vary sub¬ 
stantially across the frequency extent of a typical spectral 
line. For example, if we take the CO J = 1 — 0 transi¬ 
tion at 115 GHz, then for a spectral line with a full width 
zero intensity of 200 MHz (corresponding to an extent of 
520kms _1 ), the fractional change in the Planck function 
across the spectral line is less than 0.3% for all excitation 
temperatures above the CMB. For most sources within 
the Milky Way Galaxy, the extent of a spectral line is 
an order of magnitude smaller and the corresponding 
change in the Planck function across the transition is also 
smaller. The fractional uncertainty in the Planck func¬ 
tion increases as the frequency of the transition increases; 
however, even at 1 THz, the fractional uncertainty is still 
less than 3% for a 500km/s wide transition. As a result, 
the Jv{T ex ) — Jv(Tb g ) term is approximately constant in 
all realistic cases and the integral term in Equation [721 
reduces to just the integrated intensity of the spectral 
line: 


j\rthin 
1 ^ tot 


( 3h \ / Q rot \ ex P (m&) 
V8? T 3 Sg, 2 Ri) \gjgKgi) exp ( _ i 

y kl ex J 

x _1_ f T R dv 

(Jv{T ex ) - J„{T bg )) J f 


(80) 


N< 


thin-\- RJ 


tot 


= i- 


3 k 


\8TT 3 vSp 2 Ri 


(-2i£u) ex p(_5A 

\gjgngiJ \kT ex j 


T R T ex dv 
f ( T ex - T bg ) ■ 


(81) 


We test the validity of these approximations by compar¬ 
ing the ratio of Equation [HU] (A ) to Equation IHTI 
(N^ n+RJ ) in Figure [3j The column density calcu¬ 
lated from Equation [HU] is always less than the column 
density calculated from Equation [HU in the Rayleigh- 
Jeans limit. Darker shading corresponds to worse agree¬ 
ment between the equations. For example, the Rayleigh- 
Jeans approximation is a good approximation (less than 
5% disagreement between equations) for the low fre¬ 
quency (4.8 GHz) H 2 CO K-doublet transition for all exci¬ 
tation temperatures. The Rayleigh-Jeans approximation 
is still a good approximation for the NH 3 inversion tran¬ 
sitions (e.g.(l,l), (2,2), etc.) at 24GHz with less than 
10% disagreement between the equations for T ex > 5K. 
Rayleigh-Jeans failure is more apparent for transitions 
at millimeter and shorter wavelengths. For instance, the 
C 18 0 J = 1 — 0 line at 109.7 GHz has more than 15% 
disagreement between the optically-thin and optically- 
thin plus Rayleigh-Jeans approximations for T ex < 9K 
that rapidly increases toward lower excitation tempera¬ 
tures. It is possible to use Equation [HT] at submillimeter 
wavelengths v > 300 GHz if the excitation temperature 
is large (i.e. less than 10% disagreement at 300 GHz if 
T ex > 26 K); however, for most cases in dense regions 
of molecular clouds where the tracer is sub-thermally 
populated with low excitation temperatures, Eauation ISOI 
should be used. 


11.3. Optically Thin with Rayleigh-Jeans 
Approximation and Negligible Background 

Assuming that the temperature of the background 
source (i.e. the cosmic microwave background radiation) 
is small in comparison to the molecular excitation tem¬ 
perature ( T bg -C Tex), Equation IHU becomes: 


yy thin-\-RJ-\-noBG _ 


3k 


8n 3 vSg, 2 R., 


) ( Qrot ^ 

\gjgKgi) 


x exp 


E u 

kTe, 


gjgngi J 
T R dv 

~1~' 


(82) 


We test when it is appropriate to use this approx¬ 
imation by comparing the column density calcu¬ 
lated with Equation [HU] (N^ n ) versus Equation [571 
^jythin+ R j+n° B Gy rat j 0 0 f these two equations is 
given by N thin /N thin +R J + noBG = j^ Tex )/[J„(T ex ) - 

Ju{T bg )\. Figure [5] shows contours of this ratio as a func¬ 
tion of T ex and the transition frequency. The general 
trend is such that the thin+RJ+noBG approximation 
is poor for low frequencies and for low T ex . At first, 
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Figure 3. In the top panel we compare the molecular col¬ 
umn density calculated assuming optically-thin (A r ^ v/ "; Equa¬ 
tion [80j and optically-thin within the Rayleigh-Jeans approxima¬ 
tion ; Equation 1811) . In the bottom panel we show the 

comparison between molecular column density calculated assuming 
optically-thin and optically-thin within the Rayleigh-Jeans approx¬ 
imation and with the assumption of negligible background emission 
{N thin +R J + noBG. Equation [g2j. 

it seems counterintuitive that agreement between the 
two approximations would be better at high frequencies 
when failure of the Rayleigh-Jeans limit is more likely 
to occur (see Section H 1.2 11 : however, ignoring the back¬ 
ground term in Equation [82] compensates for Rayleigh- 
Jeans failure. For high frequency and for T ex » T bg , the 
J v {Tbg) term in the equation ratio becomes small such 
that N%™/N%i n+RJ+noBG 1. The disagreement be¬ 
tween the two equations is worse at smaller T ex . 

As a first example, we consider the NH 3 (1,1) lowest 
inversion transition at 23.7 GHz. At Tp X = 10 K, a com¬ 
mon value found in NH 3 studies (see iRosolowskv et, al.1 
I2008D . the thin+RJ+noBG approximation results in a 
systematic 30% underestimate of the column density. At 
these frequencies, the CMB cannot be ignored. This 
result indicates that a column density expression with 
background (Equation [HU] or [ED should be used to cal¬ 
culate column densities from NH 3 (1,1) observations. I 11 
contrast, the thin+RJ+noBG approximation is more ap¬ 
plicable at higher frequencies. For the HCO + J = 3—2 
transition at 267.5 GHz, Equation 1821 (thin+RJ+noBG 
approximation) has better than 10% accuracy for T ex > 
5.4 K. Equation [82] may be used for submillimeter transi¬ 
tions (v > 300 GHz) unless the level populations for the 
transition are severely sub-thermally populated ( T ex < 
5K). 


O' 


H-1-1— 


0.1 


0.6 


0.5 

-I— 


■s 

*■6 

2 


^thin ^ J\Jthin+RJ 


H-!-1--h 


Klthin / Mthin+RJ+noBG _ 

^tot / iN tot 


Better than 
1% agreement 


11.4. Not Optically Thin Approximation 

If we assume that r„ > 1, we use the radiative transfer 
equation (Equation [27]) to solve for r„: 


t v = - In 



Tr 

f[J(T ex )-J(T bg )}\ ’ 


(83) 


and insert this expression for r„ in the column density 
equation ('Equation 1321) : 


a jthick 
1 ’ tot 


3 h \ 

8ir 3 Sp 2 Ri) 



( Qrot \ ^ (zh) 

\9j9k9i) exp (*§£-)- 


Tr 

f[J(T ex )-J(T bg )\_ 


dv. 


1 

(84) 


Note that Equation [84] is quite general, applicable to all 
values of zz, r, T ex , T bg: and Tr. Furthermore, one can 
derive a simple relationship between the column density 
derived assuming tCI and those derived assuming r y( 
1 by noting that one can write N^ n (EquationJED]) using 
the radiative transfer equation ('Equation [271) as follows: 


l\rthin 
1 ^ tot 


3 h \ 

8n 3 Sp, 2 Ri) 


( Qrot \ 6X P (ife) 
\9J9K9i) ex p(^-)-l 

x f (1 — exp(—r)) dv. (85) 


Therefore, a column density calculated assuming t< 1 
is related to a column density calculated more generally 
by the following: 


N, n , = N: 


thin 


1 — exp(—r) 


( 86 ) 


The factor r/( 1 — exp(—r) ) is the “optical depth cor - 
rection factor” derived by I Goldsmith fc Langcrl ([1999D . 
Note that Equation 1551 also applies to the approximation 
cases hv -C kT ex CSection [11.21) and T bg <C T ex (111.31) . 


12 . THE CONSEQUENCES OF ASSUMING A 
CONSTANT EXCITATION TEMPERATURE 

The calculation of column density requires a determi¬ 
nation of the gas excitation temperature (Equation [20]). 
The excitation temperature is a very general concept that 
describes an energy density (in units of temperature), 
whether kinetic, radiative, rotational, vibrational, spin, 
etc. 

Since molecules are fundamentally multi-level systems, 
an excitation temperature may be defined for each tran¬ 
sition allowed via electric dipole (or magnetic dipole or 
electric quadrupole) selection rules0. Most molecules (in 
particular, dense gas tracers) in the ISM are in non-LTE 
with different gas excitation temperatures for different 
transitions. However, to simplify the problem of calcu¬ 
lating the total molecular column density from a single 


8 The condition r > 1 is often somewhat erroneously referred 
to as “optically thick”, while in fact the analysis presented in this 
section is appropriate for all conditions excluding r < 1 . 

9 You can also define an “excitation temperature” between levels 
which are not radiatively coupled - e.g. the ( 1 , 1 ) and ( 2 , 2 ) inversion 
states of NH 3 . This is usually called a rotation temperature T ro t. 
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Figure 4. The optically thin column density at 1 K km/s for tran¬ 
sitions of HCO”*”. The J = 1 — 0 (blue solid line), J = 3 — 2 (red 
dotted line), and J = 4 — 3 (green dashed line) are plotted for dif¬ 
ferent excitation temperatures in the CTEX approximation. The 
inset shows the energy level diagram for the lowest 5 energy levels 
of HCO + and the three transitions easily observable from ground- 
based observatories. 

transition, the constant excitation temperature approxi¬ 
mation (CTEX) is usually assumed where all transitions 
have the same excitation temperature 

T ex = Tg TEX \/u^l . (87) 

The optically thin column densities using Equation 1801 
are shown in Figure 0] for the three lowest energy transi¬ 
tions of the HCO + molecule that are observable from the 
ground in atmospheric windows. Note that the column 
density equation diverges for Tg TEX 2.735 K due to 
the J u {T ex ) — J„(Tb g ) term in the denominator. This 
divergence has the consequence that as the level popu¬ 
lations approach equilibrium with the CMB, it becomes 
more difficult to accurately determine the column den¬ 
sity. The plotted curves show Nthin/I or the column 
density in the optically thin limit for an integrated in¬ 
tensity of IK km/s. These curves should be multiplied 
by the observed integrated intensity and corrected for 
any optical depth effects by r/( 1 — exp(—r)). The 1 — 0 
transition is the least sensitive transition to changes in 
the excitation temperature above 3 K with a minimum at 
Tg TEX = 7.46 K. As a result, the difference between as¬ 
suming an excitation temperature between 5 K and 10 K 
results in a less than 10% uncertainty in the calculated 
column density assuming CTEX. This statement is in 
general true for the 3 mm transitions of several famous 
dense gas tracers such as HCO + , HCN, N 2 H + , and CS. If 
a lower limit to the column density is desirable, then the 
column density can be calculated at the excitation tem¬ 
perature of the minimum excitation temperature curve 
(e.g. Nthin/I > 1-03 x 10 12 cm -2 / Kkm/s for the HCO + 
J = 1 — 0 transition at Tg TEX = 7.46 K.) 

Estimating the total column density from higher en¬ 


ergy transitions is more uncertain for low excitation tem¬ 
peratures. Unfortunately, this situation does occur in 
observations of dense cores in the ISM where the dense 
gas tracer is typically bot h optically thi ck and very sub- 
thermally populated (see iShirlev ct all 12013H . For the 
HCO + J = 4 — 3 transition, the uncertainty in the total 
column density calculation between assuming an excita¬ 
tion temperature of 5K versus 10 K is a factor of 40! 
The J = 4 — 3 curve does not reach its minimum until 
Tg TEX = 43.55 K, close to the energy of the J = 4 level 
above ground ( E u /k = 42.80K). It is therefore advisable 
to observe transitions with E u /k that are comparable to 
the typical excitation temperature to minimize this un¬ 
certainty. 

Proper estimation of the gas excitation temperature 
usually requires observations of multiple transitions cou¬ 
pled with modeling of the statistical equilibrium and ra¬ 
diative transfer. In general, the excitation temperature 
is a sigmoid function that varies between the background 
radiation temperature at low density and the gas kinetic 
temperature at high density. At low gas densities, colli¬ 
sions are not as important for setting the level popula¬ 
tions as interaction with the background radiation field 
and the excitation temperature equilibrates with the ra¬ 
diation temperature. At high gas densities, collisions 
dominate in setting the level populations and the excita¬ 
tion temperature equilibrates to the gas kinetic temper¬ 
ature of the dominant collisional partner. A transition 
is said to be “thermalized” if T ex is equal to Tjf. Sub- 
thermal excitation indicates that T ex is less than T k- 
In the general case of non-LTE in the ISM, it is possi¬ 
ble that some transitions for a particular molecule may 
be thermalized (i.e. CO J = 1 — 0), while higher energy 
levels are subthermally populated (i.e. CO J = 6 — 5). 
Strong non-LTE effects are also observed in some cases 
(i.e. masers where naturally occurring population inver¬ 
sions can occur). Most molecular transitions of dense gas 
tracers observed in the ISM are subthermally populated. 

Multi-level effects can cause departures from the sig¬ 
moid shape. For example, a bottleneck due to longer 
rates of spontaneous decay, A u i ~ zz , for low J tran¬ 
sitions can cause super-thermal excitation in the J = 
1 — 0 transition where, over a narrow range of densi¬ 
ties, T ex > Tk ■ A full discussion of multi-level statistical 
equilibrium with radiative transfer is beyond the scope of 
this tutorial. There are situations where the excitation 
temperature can be directly estimated from the observa¬ 
tions. If it is known that the line is very optically thick, 
then the (1 —e _T ) term in the radiative transfer equation 
(Equation 1271) approaches 1 and the excitation temper¬ 
ature may be solved analytically with an assumption of 
the gas filling factor: 


rj-iobs _ 


hu/k 


hi 1 + 


hv/k 

Tn/f+MTcmb), 


(if T » 1) . (88) 


Molecular line ratios of well resolved hyperfine transi¬ 
tions or ratios of isotopologues may be used to determine 
the optical depth in a line by taking ratios of EauationETl 
(see Appendix [Dj, although these methods usually as¬ 
sume that the excitation temperature in each hyperfine 
line or among different isotopologue transitions are the 
same which may not always be true. 

Figure [5] shows the dependence of Equation [27] on the 
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Figure 5. The ratio of the excitation temperature at a given gas 
filling fraction compared to the excitation temperature calculated 
assuming a gas filling fraction / = 1. The four curves correspond to 
spectral lines with Tr = 0.1,1,10,100 K. The curves are calculated 
for a transition at 3 mm. The curves are a weak function of the 
transition frequency. The horizontal dashed line correspond to an 
error in T ex of 10% compared to the / = 1 assumption. 

filling fraction by plotting the ratio of the excitation tem¬ 
perature at a filling fraction, /, divided by the excitation 
temperature with a filling fraction of / = 1 (the usual as¬ 
sumption when filling fraction is unknown). For intense 
spectral lines Tr > IK, calculation of T ex from Equa¬ 
tion EH is very sensitive to the gas filing fraction. For a 
line with Tr = IK, T ex is 28% larger for / = 0.5. T ex {f) 
quickly increases toward lower gas filling fractions. The 
dependence on gas filling fraction is steeper for brighter 
lines and flatter for weaker lines. The curves plotted in 
Figure [5] are not a strong function of the wavelength for 
transitions with A > 1 cm. 

If the line is optically thin, then it is not possible to 
determine the excitation temperature from a single tran¬ 
sition as the observed line temperature, Tr, is degenerate 
with excitation temperature, optical depth, and filling 
fraction. However, in this case, the excitation tempera¬ 
ture can still be estimated from multiple transitions using 
a rotation diagram. A detailed discussion of the limits of 
the CTEX approximation for HCO + observations may 
be found in Shirley et al.l (120131) . 

13. NON-LTE APPROACHES TO THE 
DERIVATION OF MOLECULAR COLUMN 
DENSITY 

Although beyond the scope of this tutorial, we should 
mention that there are a number of non-LTE approaches 
to the solution of the coupled multi-level statistical equi¬ 
librium and radiative transfer equations. These model- 
based solutions normally rely on the measurements of a 
number of molecular spectral lines from a given species 
to constrain the physical conditions within the dense 


gas environment under study. If one is interested in 
the global physical conditions within a molecular cloud, 
the radiative transfer solution for these models is sim¬ 
plified greatly by the introduction of a geometrically- 
averaged escape probability (usually denoted by /?). A 
number of escape probability expressions which depend 
upon the assumed molecular cloud geometry and global 
velocity field can be used in the se global cloud mod¬ 
els (see Ivan der Tak et all (120071) for an excellent dis¬ 
cussion). One of the most popular of these non-LTE 
escape probability models incorpor ates the Large Veloc¬ 
ity Gradient (LVG) approximation ( Sobolevlll960 ) to the 
solution of the radiative transfer and statistical equilib¬ 
rium equations. At the heart of the LVG approxima¬ 
tion is the assumption that photons which are emitted 
within the modelled molecular cloud escape quickly from 
the cloud due to the large velocity gradient (often as¬ 
su med to be a result of cloud collap se). As described 
in iMangum fe Woottenl (ll993l) and Ivan der Tak et all 
ms: in the LVG approximation the escape probability 
is given by /3 = —|— within a uniform collapsing spher¬ 
ical geometry. The molecular cloud is assumed to have 
uniform volume density, kinetic temperature, and molec¬ 
ular column density. Inputs to an LVG model include a 
specification of the molecular spectrum and collisional 
excitation rates involving the primary collision partner 
in the molecular cloud (usually assumed to be H 2 ) and 
the molecular species under study. Although the LVG 
approximation entails the use of a rather simplified uni¬ 
form physical model, it often produces plausible physical 
conditions for environments where the transitions used 
to constrain the model are not too optically thick. 

Other non-LTE approaches i nclude the generalized es - 
cape probability code RADEX (Ivan der Tak et al . 20071) 
microturbulent velocity models dLeung fc Liszt! 19761) 
and Monte Carlo approaches (jBrinch fc Hogerhciidd 
120101b In the microturbulent model solution the molec¬ 
ular cloud is assumed to have a power-law distribution 
in volume density and kinetic temperature with a line 
profile that is produced due to a combination of ther¬ 
mal and turbulent motions. The radiative transfer equa¬ 
tion is iteratively solved using a quasi-diffusion method. 
Monte Carlo models are generalized solutions to the ra¬ 
diative transfer in molecular clouds which incorporate 
up to three dimensions of arbitrarily complex structure. 
This complex structure is organized into a grid of cells 
with uniform volume density, kinetic temperature, and 
molecular abundance. By ray tracing along random 
paths through this grid the radiative transfer equation 
is iteratively solved. 

Even though non-LTE approaches result is a more 
accurate representation of the physical conditions in a 
molecular cloud, they must be constrained by multi¬ 
ple molecular spectral line measurements from a given 
molecular species. Non-LTE models also require knowl¬ 
edge of the collisional excitation rates for the molecu¬ 
lar species under study. Collisional excitation rates for 
a wide variety of molecules have been calculated (see, 
for example , the Leiden A t omic and Molecular Database 
(LAMDA); ISchdier et al.l (I2005IH . In general, though, 
limitations to their accuracy and kinetic temperature 
range can result in added uncertainty in non-LTE mod¬ 
elling of radiative transfer. Due to these limitations, LTE 
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calculations of the molecular column density are often 
satisfactory with limited measurements. 


14. BEAM-AVERAGED VERSUS 
SOURCE-AVERAGED COLUMN DENSITY 

An aspect of the calculation of molecular column den¬ 
sity from a spectral line measurement that is often 
overlooked is the difference between beam- and source- 
averaged measurements. Beam-averaged molecular col¬ 
umn densities are derived from spectral line measure¬ 
ments with a defined spatial resolution element (a radio 
telescope “beam” 9b ea m)■ Without an understanding of 
the actual spatial extent of the measured emission (i.e. 
Obeam > 0source), one normally assumes that it is simply 
averaged over the beam of the measurement, or “beam- 
averaged” . For measurements which resolve the spatial 
extent of emission (i.e. 9b ea m < 0source), spectral line 
measurements measured over a source, usually consid¬ 
ered “extended”, result in “source-averaged” molecular 
column densities. 


15. MOLECULAR COLUMN DENSITY 
CALCULATION EXAMPLES 

In the following we describe in detail some illustra¬ 
tive calculations of the molecular column density. These 
example molecules were chosen for their commonality 
(C 18 0 is an abundant, though often optically-thin, iso- 
topomer of CO), relatively simple hyperfine structure 
(C 17 0), somewhat more complex hyperfine structure 
(N 2 H+), or general usage as kinetic temperature (NH 3 ) 
or volume density (H 2 CO) probes. In these examples 
we adopt the more conventional notation for a generic 
(emission or absorption) transition, i.e. “ J = 1 — 0”, in¬ 
stead of the more explicit (and somewhat less convenient) 
“ J = 1 —► 0” notation. We also address a common point 
of confusion for students learning this material and its re¬ 
liance on CGS (Centimeter-Gram-Second) units. For the 
calculations presented below dipole moments are quoted 
in the CGS unit “Debye”, which is equal to 10 “ 18 esu cm 
(electrostatic units times centimeter). A good source for 
molecular dipole moments ( ft ) and rotational constants 
(i.e. A 0 , B 0l Cq) is the Jet Propulsion Laboratory (JPL) 
Mol ecular Spectroscopy database and spectral line cata¬ 
log ([Pickett ct aLlllblMFl 


15.1. C 18 0 

To derive the column density for C 18 0 from a measure¬ 
ment of its J = 1 — 0 transition we use the general equa¬ 
tion for molecular column density (Equation [32]) with the 


10 Available online at http://spec.jpl.nasa.gov (use the “cata- 
log directory with links” to access dipole moment and molecular 
constants information). 


following properties of the C 18 0 J = 1 — 0 transition: 

S — ^ u 

2 J u + 1 

ft, = 0.11079 Debye = 1.1079 x 10~ 19 esu cm 
B 0 = 54891.420 MHz 
gj = 2 J u + 1 

gpc = 1 (for linear molecules) 
gi = 1 (for linear molecules) 
kT 1 

Qrot — tt 7 + n (Equation O 
liB 3 

~ 0.38 (T + 0.88) 

E u = 5.27 K 
v = 109.782182 GHz 

which leads to: 

3 h 


N tot (C ls O) = 


8n 3 fj, 2 J u Ri 
exp 


kT ex 1 
~hB + 3 I exp 


Eu 

kT., 


hu 


kT., 


- 1 


J r v dv. (89) 


Assuming that the emission is optically thin (t v -C 1; 
Equation [SOll. Equation IM1 becomes: 


N to t{C l «0) = 


2.48 x 10 14 (T ex + 0.88) exp (f^) 


exp 


(«) 


- 1 


J Tndv(km/s) 
lf{JATe X )-MT bg ))\ 


cm 


-2 


(90) 


If we are using integrated fluxes (S„Av) instead of in¬ 
tegrated brightness temperatures, we use Equation IB2I 
(see Appendix iBl) . with the assumption that t<1): 


NtotiC^O) = 


3c 2 


( Qrot j 


1 6Tr 3 n s Sfj, 2 iy 3 \gjgKgiJ 

En, 


x exp 


S u Av 


kT 

4.77 x 10 15 (Tex + 0. 


exp 


(«?) 


@maj (■ asec)9 rn i n (asec ) 


exp 


(w) 


f S v (Jy)dv(km/s) 
[f(MT ex )-MT bg ))\ 


cm 


-2 


- 1 
(91) 


where we have assumed a Gaussian source (f2 s = 
1.133 dmajOmin) with 9 maj and 9 min the major and mi¬ 
nor source diameters (in arcseconds). 


15.2. C 17 0 

C 17 0 is a linear molecule with hyperfine structure due 
to interaction with the electric quadrupole moment of 
the 17 0 (I = |) nucleus. Using the selection rule: 

F = J + I, J + I — 1, J + I — 2, — I\, 

we find that each J-level is split into the hyperfine levels 
indicated in Table |4] (for the first five J-levels). Since 
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Table 4 

Allowed C 17 0 Hyperfine Energy Levels 


J 

Number of Energy Levels 

Allowed F 

0 

i 

5 

? 5 3 

1 

3 

2 

5 

§’ ?’ § 3 1 

3 

6 

?f 2 9 2 7 2 S 2 3 1 
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7 

S'?l 2 9 2 7 2 5 2 3 1 

2’ 2 ’ 2 ’ 2 ’ 2 ’ 2 ’ 2 


J 
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c 17 o 
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5/2 


3/2 
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1 


< 
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3/2 
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Figure 6. Electric quadrupole hyperfine energy level structure for 
the J=0, 1, and 2 levels of C 17 0. Note that the 3 (J = 1 — 0) and 
9 (J = 2 — 1) allowed transitions are marked with arrows ordered 
by increasing frequency from left to right. 


Av (km/s) 
0 



Av (km/s) 

-10 1 



Figure 7. Synthetic spectra for the C 17 0 J = 1 — 0 (top) and 
2 — 1 (bottom) transitions. Horizontal axes are offset velocity (top) 
and frequency (bottom) relative to 112359285.0 and 224714386.0 
kHz, respectively. Overlain as a dashed line is a synthetic 100 kHz 
gaussian linewidth source spectrum. 


the selection rules for the single-spin coupling case is, 
A F = 0, ±1, and A J = ±1, there are 3, 9, and 14 allowed 
hyperfine transitions for the J = 1 — 0, J = 2 — 1, and 
J = 3 — 2 transitions, respectively. Figure [6] shows the 
energy level structure for the J = 1 — 0 and J = 2 — 1 
transitions. 

We can calculate the relative hyperfine intensities ( Ri ) 
for the J = 1 — 0 and J = 2 — 1 transitions using 
the formalism derived in Section [TOj Using the formal¬ 
ism described in Section [10] we can derive the relevant 
Ri for the electric quadrupole hyperfine coupling cases 
(. Ri(F,J),I = |; Table O. Figure [7] shows the synthetic 
spectra for the C 17 0 J = 1 — 0 and J = 2 — 1 transitions. 

To derive the column density for C 17 0 from a measure¬ 
ment of its J = 1 — 0 transition we use the general equa¬ 
tion for molecular column density (Equation 1321) with 
the following properties of the C^O J = 1 — 0F=| — | 


transition: 


2 J u + 1 

fi = 0.11032 Debye = 1.1032 x 10~ 19 esu cm 
S 0 = 56179.9911 MHz 
gj = 2 J u + 1 

gx = 1 (for linear molecules) 
gi = 1 (for linear molecules) 
kT 1 

Qrot (Equation [52]) 

hB 3 

~ 0.37 (T + 0.90) 

S„ = 5.39 K 














































How to Calculate Molecular Column Density 


17 


Table 5 

Hyperfine Intensities 3, for C 17 0 J = 1 — 0 and J = 2 — 1 


F' -> F 

J' -4- J 

a 

b 

c 

(2F' + 1)(2F+1) 
(27+1) 

Au h (kHz) 

6j 

R i(F,J) 

(3 5n 

V 2 ’ 2 ' 

(1,0) 

5 

2 

0 

5 

2 

4 

-501 

1 

3^2 
_ 3+2 

2 

9 

(1 5\ 
v 2 ’ 2 ' 

(1,0) 

5 

2 

1 

7 

2 

8 

-293 

4 

9 

(- ~) 

V2 ’ 2' 

(1,0) 

5 

2 

5 

2 

i 

6 

+724 

1 

3v/2 

3 

9 

(- ~) 

V 2 ’ 2 ' 

(2,1) 

b 

2 

i 

t) 

2 

4 

-867 

i 

10 _ 

1 

25 

(- -) 

v 2 ’ 2 ' 

(2,1) 

5 

2 

5 

2 

2 

6 

-323 

_ jV2_ 

15^7 

64 

525 

(1 5) 

^2 ’ 2 ' 

(2,1) 

5 

2 

2 

7 

2 

8 

-213 

V3 

2y/35 

6 

35 

1) 

V2 ’ 2 / 

(2,1) 

5 

2 

2 

9 

2 

40 

3 

-169 

1 

2 VT 0 

1 

3 

(- ~) 

v 2 » 2 / 

(2,1) 

5 

2 

1 

3 

2 

4 

3 

-154 

1 

2y/5 

1 

15 

( — — ) 

(2,1) 

5 

2 

3 

2 

2 

8 

3 

+358 

+7 

10+2 

7 

75 

(- ~) 

v 2 ’ 2 ' 

(2,1) 

5 

2 

i 

7 

2 

8 

+694 

_ 1 
6\/l4 

1 

63 

i 1 !) 

^2 ’ 2 ' 

(2,1) 

5 

2 

7 

2 

2 

32 

3 

+804 

1 

4V7 

2 

21 

(5 3 \ 

V2 ’ 2^ 

(2,1) 

5 

2 

2 

5 

2 

4 

+902 

s/7 

15+2 

14 

225 


3 The sum of the relative intensities J+ ft., = 1.0 for each A J = 1 transition. 
b Frequency offsets in kHz relative to 112359.285 and 224714.368 MHz for J = 1 — 0 
and J = 2 — 1. respectively (from the spectroscopic constants in IKlappcr et all (120031) ). 


which leads to: 

N tot (C ir O) = 


3 h 


kT P , 
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8 t T 3 ii 2 J u Ri V hB 3 

n -1 


exp 


E u 

kT e , 


exp 


hu 

kT„, 


- 1 


J r u dv. (92) 


Assuming that the emission is optically thin (t„ -C 1; 
Equation 1501). Equation [Ml becomes: 


N to t(C 17 0) 


2.44 x 10 14 (T ex + 0.90) 

R-i 

f Tfidv(km/s) 

< [f(MT eX )~ MTbg)) 


ex p (fcfc) 

ex p (fcfe) - 1 

■ cm -2 , (93) 


where v is the frequency of the hyperfine transition 
used. For example, if the F=| — | hyperfine was cho¬ 
sen for this calculation, R, = | (See Table [5]) and 
z/ = 112359.285 + 0.724 MHz = 112.360009GHz. Equa¬ 
tion [53] then becomes: 


Ntot(C 17 0) 


= 7.32 x 10 14 (T ex + 0.90) 



/ Tndv(km./s) 

J (MT ex ) - MT bg ))_ 


cm 


(94) 


15.3. N 2 H+ 

N 2 H+ is a multiple spin coupling molecule due to the 
interaction between its spin and the quadrupole moments 
of the two nitrogen nuclei. For a nice detailed descrip- 
tio n of the hyperfine le v els of the J = 1 — 0 transition 
see iShirlev et al.l (120051) : iPaeani et al.1 (|2009f) . Since the 
outer N nucleus has a much larger coupling strength than 


Table 6 

Outer Nitrogen ( F\) Hyperfine Intensities for N 2 H+ J = 1 — 0 


+(-+ +i a 

J' -+ J 3 

a 

b 

c 

(2F(+l)(2F 1 + l) 
(2+v+l) 

6 j 

Ri(+i, J) 

( 0 , 1 ) 

(i,i) 

( 1 , 0 ) 

( 1 , 0 ) 

i 

i 

0 

1 

1 

1 

1 

3 

1 

5 i 

3 

1 

! 

( 2 , 1 ) 

( 1 , 0 ) 

i 

1 

2 

5 

1 

3 

9 


3 In = 1- 


the inner N nucleus, the hyperfine structure can be deter¬ 
mined by a sequential application of the spin coupling: 

Fi = J + In , 

F = F\ + In- 

When the coupling from both N nuclei is considered: 

• The J = 0 level is split into 3 energy levels, 

• The J = 1 level is split into 7 energy levels, 

• The J = 2 and higher levels are split into nine 
energy levels. 

Since the selection rules for the single-spin coupling case 
apply, AFi = 0, ±1, A F = 0, ±1, 0 0, and A J = ±1, 

there are 15, 34, and 38 hyperfine transitions for the J = 
1 — 0, J = 2 — 1, and J = 3 — 2 transitions, respectively. 
Figure[5]shows the energy level structure for the J = 1 — 0 
transition. 

To illustrate the hyperfine intensity calculation for 
N 2 H+, we derive the relative intensities for the J = 1 — 0 
transition. Relative intensities, derived using the formal¬ 
ism presented in Section [Tol and Equation [fT] are listed 
in Tables [ 6 ] [7] and[S] Figure [5] shows the synthetic spec¬ 
trum for the N 2 H + J = 1 — 0 transition. 

To derive the column density for N 2 H + we start with 
the general equation for the total molecular column den- 
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Table 7 

Inner Nitrogen (F) Hyperfine Intensities for N 2 H+ J = 1 — 0 


F' F a 

F[ -> Fi a 

a b 

c 

(2F' + 1)(2F+1) 
(2/w+l) 

6 j 

RdSFi) 

(1.0) 

(1.1) 

(0,1) 

(0,1) 

i 

i 

0 

1 

1 

1 

1 

3 

1 

3 i 

,3 

1 

| 

(b2) 

(0,1) 

i 

1 

2 

5 

1 

? 

3 i 

3 

I 

? 

! 

12 

(0,1) 

(1,0) 

(1,1) 

(1.1) 

i 

i 

1 

1 

1 

1 

1 

1 

(1,1) 

(1,1) 

i 

1 

1 

3 

1 

6 

(1,2) 

(1,1) 

i 

1 

2 

5 

1 

6 

f 

¥ 

12 

(2,1) 

(2,2) 

(1.1) 

(1,1) 

i 

i 

1 

2 

2 

1 

5 

25 

3 

1 

6 i 

2v/6 

(1,0) 

(2,1) 

i 

2 

1 

i 

1 

? 

6 

1 

! 

¥ 

1J0 

4 

(1,1) 

(2,1) 

i 

1 

2 

3 

(1,2) 

(2,1) 

i 

1 

2 

5 

1 

30 

(2,1) 

(2,1) 

i 

2 

2 

5 

!. 

2^ 

5 

(2,2) 

(3,2) 

(2,1) 

(2,1) 

i 

i 

2 

2 

2 

3 

25 

A 

3 

1 

¥ 

15 


a In = 1- 


Table 8 

Hyperfine Intensities 3, for N 2 H+ J=1 — 0 


F' ->• F b 

F[ -> Fi b 

J' -> J 

RiCFi, JJRiOF.Fi) 

Au h (kHz) 

Ri (obs) c 

(0,1) 

(1,1) 

(1,0) 

1 

27 

-2155.7 

l 

27 

(2,2) 

(1,1) 

(1,0) 

% 

-1859.9 

z 

27 

(2,1) 

(1,1) 

(1,0) 

5 

108 


(1,2) 

(1,1) 

(1,0) 

5 

108 

-1723.4 

I 

9 

(1,1) 

(1,1) 

(1,0) 

1 

36 



(1,0) 

(1,1) 

(1,0) 

1 

27 



(2,1) 

(2,1) 

(1,0) 

% 

3$ 

-297.1 

5 

27 

(2,2) 

(2,1) 

(1,0) 

5 

108 


(3,2) 

(2,1) 

(1,0) 

27 

+0.0 

7 

27 

(1,1) 

(2,1) 

(1,0) 

Z 

108 

+189.9 

1 

9 

(1,2) 

(2,1) 

(1,0) 

1 

324 



(1,0) 

(2,1) 

(1,0) 

5 

V 



(1,2) 

(0,1) 

(1,0) 

5 

8,1 

+2488.3 

1 

9 

(1,1) 

(0,1) 

(1,0) 

1 

2.7 



(1,0) 

(0,1) 

(1,0) 

1 

81 




a The sum of the relative intensities Ri = 1.0. 

b Frequency offset in kHz relative to 93173.7767 MHz dCaselli et al.|[T9950 . 
c Since the J=0 level has no hyperfine splitting, only the sum of all transitions 
into the J=0 is observed. 


sity ('Equation 13^1) with: 

S = J u (see Section 19.11) 

4“ 1 

/i = 3.40 Debye = 3.40 x 10” 18 esu cm 
S 0 = 46586.88 MHz 

Ri = (see Section [lU] or, for J=l-0, see Table 0 


9J 

9K 

9i 


Qrot 


2 J u + 1 

1 (for linear molecules) 
1 (for linear molecules) 
kT 1 

— + - (Equation [52j 

tit) 6 

0.45 (T + 0.75) 


K = 4.4716 K 


which leads to: 




exp 


hv 

kT e , 


E u 

kT P . 


- 1 


x 


-l 


TydV. 


(95) 
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Fi 


I^H + 


Outer Nitrogen Coupling 


Inner Nitrogen Coupling 


Figure 8. Energy level structure for the J — 1 — 0 transition 
of N 2 H+. Note that of the 15 hyperfine split levels only 7 are 
observed due to the fact that there is no hyperfine splitting of the 
J=0 level, but it is degenerate in energy. Grouping of the indicated 
transitions show the 7 observed transitions. Transitions are ordered 
by increasing frequency from left to right. 



hkV (kHz) 


2000 


Figure 9. Synthetic spectra for the N 2 H + J = 1 — 0 transition. 
Horizontal axes are offset velocity (top) and frequency (bottom) 
relative to 93173.7767 MHz. Overlain as a dashed line is a synthetic 
100 kHz (Av = 0.3218 kms -1 ) gaussian linewidth source spectrum. 


Assuming optically thin emission (Equation Ej), we find 
that Equation 1551 becomcJ n l: 


N tot (N 2 H+) = 


3.10 x 10 11 (T ex + 0.75) 


exp 


(fcT ex ) 


f Tndv(km/s ) 

f (MTex) - MTbg)) J 


exp 


( hu \ 

\ kTex J 


~ 1 


cm 


(96) 


11 Alternatively, one can derive r using the ratio of two hyperfine 
transition intensities as described in Appendix [D] 


where v is the frequency of the hyperfine transition used. 
For example, if the F=(2,l), J=(1,0) hyperfine was cho- 
sen for this calculation, i?, ; = jp (See Table 0) and 
v = 93.1737767GHz. Equation l96l then becomes: 

N tot (N 2 H+{ 3, 2; 2,1)) = 1.20 x 10 12 (T ex + 0.75) 

^ ’ fT R dv(km/s) 
exp (Ml) - 1 X If (MTex)~ MT bg )) 

where the transition designation is in the (F'jF'^FjFi) 
format. Note that this transition is blended with the 
nearby F' = 1 and 2 (F) .F,) = (2,1) hyperfine transitions 
under many molecular cloud physical conditions. 

15.4. NH 3 

Ammonia (NH 3 ) is a symmetric top molecule with 
three opposing identical H (spin=l) nuclei. Quantum 
mechanical tunneling of the N nucleus through the po¬ 
tential energy plane formed by the H nuclei leads to in¬ 
version splitting of each NH 3 energy level. On top of this 
inversion splitting the energy levels are split due to two 
hyperfine interactions: 

J—Iat:: Coupling between the quadrupole moment of the 
N nucleus (Ijv = 1) and the electric field of the 
H atoms, which splits each energy level into three 
hyperfine states. For this interaction the angular 
momentum vectors are defined as follows: Fi = 
J + In- 

Fi— Ih" Coupling between the magnetic dipole of the 
three H nuclei (I# = I for Ky 0 or 3n (para), 
Iff = | for K=3n (ortho)) with the weak current 
generated by the rotation of the molecule. For this 
interaction the angular momentum vectors are de¬ 
fined as follows: F = F\ + Iff. 

Weaker N-H spin-spin and H-H spin-spin interactions 
also exist, but only represent small perturbations of the 
existing hyperfine energy levels. 

Figure [TU] shows all NH 3 energy levels below 1600 K, 
while Table [9] lists the level energies for J < 6 . Figure fill 
shows the rotational energy level diagram for the first 
four J-levels of NH 3 , while Figure El shows the inver¬ 
sion and hyperfine level structure for the ( 1 , 1 ) and ( 2 , 2 ) 
transitions. 

We can calculate the relative hyperfine intensities () 
for the (1,1), (2,2), (3,3), and (4,4) transitions using the 
formalism derived in Section 1101 Using the formalism 
presented in Section [TU] we can derive the relevant Ri for 
the quadrupole hyperfine (F^(F 3 ,T), Iv — 1: Tables El 
El El andfTTl) and magnetic hyperfine (7?,;(F, Fj), I H = 
I (para, K ^ 0 or 3n) or | (ortho, K = 3n); Tables fl5l 
El El and E| coupling cases. The resultant total hy¬ 
perfine intensities are derived by multiplying the mag¬ 
netic hyperfine intensities (i?,(F, F 3 )) by the associated 
quadrupole hyperfine intensities (Ri(Fi, J)). Tables [TUI 
and El show the results of this calculation for the NH 3 
(1,1) and ( 2 , 2 ) transitions, respectively. F igure [TUI shows 
the synthetic spectra for the NH 3 (1,1) and (2,2) transi¬ 
tions. 


cm 2 . (97) 
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Figure 10. Rotational energy level diagram for NH 3 . All levels with energy < 1600 K are shown. 


Table 9 



NH 3 Level Energies 3 

,b 

O 

- 

Level 

Energy (K) 

Level 

Energy (K) 

GO 


( 0 , 0 ,a) 

1.14 




- 

( 1 , 1 ,s) 

23.21 

( 1 , 1 ,a) 

24.35 



( 1 , 0 ,s) 

28.64 





( 2 , 2 ,s) 

64.20 

( 2 , 2 ,a) 

65.34 

CO 

— 

( 2 , 1 , 8 ) 

80.47 

( 2 , 1 ,a) 

81.58 


- 

( 2 , 0 ,a) 

86.99 





(3,3,s) 

122.97 

(3,3,a) 

124.11 

ST 


(3,2,s) 

150.06 

(3,2,a) 

151.16 


" 

(3,1,b) 

166.29 

(3,1,a) 

167.36 

c* 9 

_ 

(3,0,s) 

171.70 



s 

H 


(4,4,s) 

199.51 

(4,4,a) 

200.66 


(4,3,b) 

237.40 

(4,3,a) 

238.48 


- 

(4,2,b) 

264.41 

(4,2,a) 

265.45 


_ 

(4,1,8) 

280.58 

(4,1,a) 

281.60 

O 


(4,0,a) 

286.98 



C\2 


(5,5,s) 

293.82 

(5,5,a) 

295.00 


- 

(5,4,s) 

342.49 

(5,4,a) 

343.58 


_ 

(5,3,s) 

380.23 

(5,3,a) 

381.25 



(5,2,s) 

407.12 

(5,2,s) 

408.10 



(5,l,s) 

(5,0,s) 

( 6 , 6 ,s) 

423.23 

428.60 

(5,1,a) 

424.18 

O 


405.91 

( 6 , 6 ,a) 

407.12 



(6,3,s) 

551.30 

( 6 ,3,a) 

552.25 



( 6 , 0 ,a) 

600.30 






- 20 - 


NFL 


10 
10 " 


x 

o 

CO 

05 

cvi 

10 


21 


23.099 GHz 


10 

23.722 GHz = 
10 


ir 


23.694 GHz 


0 « 


ortho 

K = 0 


para 

K = 1 


para 

K = 2 


a Listed in level energy order per J and 
inversion-paired as appropriate. 
b See lPovnter fe Kakaii <119751) for lower-state 
energy calculations. 


Figure 11. Ammonia energy levels below 100 K. The notation 
Jx is used where d= refers to the symmetry of the energy level. 
Allowed electric dipole transitions are indicated with frequencies 
of the transitions. The total statistical weight of each energy level 
is indicated in red to the left of each energy level. 


For illustration we can derive the column density equa¬ 
tion for a para-NH 3 (Ky^O or 3n) inversion (A K = 0) 
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NH, (1,1) 


F, 



< 


< 



Magnetic Dipole F = Fi + Ig 


F 

3/2 

1/2 

3/2 

5/2 


1/2 


1/2 

3/2 

3/2 

5/2 


1/2 


NH, (2,2) 





< 


F 

3/2 

5/2 


< 


5/2 

7/2 


< 


1/2 

3/2 



3/2 

5/2 

5/2 

7/2 


< 


tt 


Magnetic Dipole F = Fi + Ih 


1/2 

3/2 


Figure 12. Inversion and hyperfine energy level structure for the (1,1) (top) and (2,2) (bottom) transitions of NH 3 . Note th at the 18 (1,1) 
and 2 4 (2,2) allowed transitions are marked with arrows ordered by increasing frequency from left to right. Adapted from [Ho &; Townes! 
(fl983h . 
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Av (km/s) 
0 


20 
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0 
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Figure 13. Synthetic spectra for the NH 3 (1,1) (top) and 
(2,2) (bottom) transitions. Horizontal axes are offset veloc¬ 
ity (top) and frequency (bottom) relative to 23694.495487 and 
23722.633335 MHz, respectively. Overlain as a dashed line is a syn¬ 
thetic 100kHz (Av = 1.265 and 1.264kms —1 , respectively) gaus- 
sian linewidth source spectrum. 


transition. For para-NH 3 inversion transitions: 


Ju(Ju + 1 ) 

/i = 1.468 Debye = 1.468 x 10 -18 esu cm 
Ri = (see Section [TU] or, for ( 1 , 1 ) and ( 2 , 2 ), see 
Tables [T9] and IUU1) 
gj = 2 J u + 1 
g K = 2 for 0 

9i = \ for 3n 

Doing the calculation in two stages, we can first derive 
the following equation for the molecular column den¬ 
sity in NH 3 as derived from a measurement of a (J,K) 
inversion (AJ = 0 ,A K = 0 ) transition assumi ng opti¬ 
cally thin emission (see iRosolowskv et all (120081) for the 
methodology to fit hyperfine transitions in the optically 


thick limit) using Equation [ 


N tot (NH 3 ) = 


3 h J u (Ju + 1) Q rot 

8TT 3 g 2 Ri K 2 gjgKgi 


exp 1 


kT ex J 


exp 


( hv ^ 
\ kT ex J 


- 1 


f T R dv 


f (MT ex ) - MT bg )) 


(98) 


Further assuming that we are measuring para-NH 3 (K^0 
or 3n), Equation l98l becomes: 


N tot (NH 3 ) ~ 


exp 


\ kT ex J 


7.44 x 10 12 J„(J U + l)Q ro t 

RiK 2 (2J u + 1) 

f T R dv(km/s) 


exp 


( hv ^ 
y kTex J 


~ 1 


f (MTex) - MT bg )) J 


cm 


(99) 


15.5. H 2 CO 

Formaldehyde (H 2 CO) is a slightly asymmetric ro¬ 
tor molecule with k ~ —0.961 (Equation l45l) . which 
means that H 2 CO is nearly a prolate symmetric rotor. 
The slight asymmetry in H 2 CO results in limiting pro¬ 
late (quantum number K_i) and oblate (quantum num¬ 
ber K + i) symmetric rotor energy levels that are closely 
spaced in energy, a feature commonly referred to as “K- 
doublet splitting”. Figure [14] shows the energy level dia¬ 
gram for H 2 CO including all energy levels E < 300 K. In 
addition to the asymmetric rotor energy level structure 
H 2 CO possesses spin-rotation and spin-spin hyperfine en¬ 
ergy level structure. Magnetic dipole interaction between 
the H nuclei (1^ = 0/1 for para-/ortho-H 2 CO) and ro¬ 
tational motion of the molecule result in spin-rotation 
hyperfine energy level splitting. Since the two H nu¬ 
clei have equal coupling, the two H nuclei spins couple 
with each other first and the hyperfine coupling scheme 
is Ih = h + h and F = J + I R ■ Since the interchange 
of the two 1 H atoms (/1 = J 2 = 1/2) obey Fermi-Dirac 
statistics (total wave function is anti-symmetric to this 
interchange) and alternating K a ladders in H 2 CO have a 
different rotational wavefunction symmetry, only ortho- 
112 CO levels ( Ih = 1) with K a = odd integers have hy¬ 
perfine splitting. For the ortho-H 2 CO lio — In transi¬ 
tion, the frequency offsets of these hyperfine transitions 
are A v < 18.5 kHz. The weaker spin-spin interactions 
between the nuclei are generally not considered. 

Table [TU] lists the frequencies and relative intensities 
for the spin-rotation hyperfine transitions of the H 2 CO 
lio — I 11 , 2n — 2i2, and 3 i 2 — 3i3 transitions. Fig- 
urefTKlshows the synthetic spectra for the H 2 CO lio—In, 
2n — 2i2, and 3 i 2 — 3i3 transitions. Note that the hy¬ 
perfine intensities are exactly equal to those calculated 
for the spin-rotation hyperfine components of NH 3 (see 
Section [0. 

For illustration we can derive the column density equa¬ 
tion for a ortho-H 2 CO [K -1 odd) K-doublet (AiF_i = 0) 
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Figure 14. Rotational energy level diagram for H 2 CO including all energy levels with E < 300 K. 


Table 10 

Fi-I# Hyperfine Frequencies and Intensities for H 2 CO J=l, 2, and 3 K-Doublet 

Transitions 


F' -> F a 

J' -> J a 

A HF h (kHz) 

a 

b 

c 

(2F'+1)(2F!+1) 

(2/jf+l) 

6j 

Ri (+ 1 ,J) 

(1,0) 

(1,1) 

-18.53 

1 

1 

1 

1 

1 

? 

_ ? 

2+5 

1 

? 

§ 

12 

(0,1) 

(1,1) 

-1.34 

1 

1 

1 

1 

(2,2) 

(1,1) 

-0.35 

1 

2 

1 

25 

3 

(2,1) 

(1,1) 

+4.05 

1 

1 

2 

5 

1 

6 

5 

ae 

(1,2) 

(1,1) 

+6.48 

1 

1 

2 

5 

1 

6 

7> 

3,6 

(1,1) 

(1,1) 

+11.08 

1 

1 

1 

3 

1 

6 

1 

¥ 

20 

(1,1) 

(2,2) 

-20.73 

1 

1 

2 

3 

1 

2 Vp 

(1,2) 

(2,2) 

-8.5 

1 

2 

2 

5 

1 

\° 

10 

1 

2 i° 

20 

(2,1) 

(2,2) 

-0.71 

1 

2 

2 

5 

(3,3) 

(2,2) 

+0.71 

1 

3 

2 

49 

3 

_ 2+2 
3+35 

392 

945 

(3,2) 

(2,2) 

+1.42 

1 

2 

3 

35 

4 

2^3 

& 

3 

1 

15 

7 

1^5 

T 35 

(2,3) 

(2,2) 

+9.76 

1 

2 

3 

1 

¥ 

6 i 

P 

28 

(2,2) 

(2,3) 

(2,2) 

(3,3) 

+10.12 

1 

1 

2 

3 

2 

3 

25 

¥ 8 

189 

(4,3) 

(3,3) 


1 

3 

4 

21 

3 

112 

(4,4) 

(3,3) 

+0.00 

1 

4 

3 

27 

+5 

4+2T 

84 

45 

112 

(3,3) 

(3,3) 

-10.4 

1 

3 

3 

49 

3 

121 

432 

(2,2) 

(3,3) 

+23.0 

1 

2 

3 

25 

3 

- 2+2 

3+35 

40 

189 

(3,4) 

(3,3) 


1 

3 

4 

21 

1 

28 

3 

H 2 

189 

(3,2) 

(3,3) 


1 

3 

3 

35 

3 

1 

21 


a I H = 1 (ortho-H 2 CO). 

b Frequency offset in kHz relative to 4829.6596 MHz for lio — I 11 , 14488.65 MHz for 2n — 
2i2, and 28974.85 MHz for 3 i2 - 3 i 3 . 
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Figure 15. Synthetic spectra for the H 2 CO lio — In, 2n — 2 i 2 , 
and 3 i 2 — 3 i 3 transitions. Horizontal axes are offset velocity (top) 
and frequency (bottom) relative to 4829.6596 MHz for lio — I 11 , 
14488.650 MHz for 2 n - 2 i 2 , and 28974.85 MHz for 3 i2 - 3i 3 , re¬ 
spectively. For the 3 i 2 — 3 i 3 transition only the A F = 0 hyperfine 
transitions are shown. Overlain as a dashed line is a synthetic 
10 kHz (Av = 0.621, 0.207, and 0.103 km s -1 , respectively) gaus- 
sian linewidth source spectrum. 


transition. For ortho-H^CO transitions: 
K 2 


S = 


Ju{Ju + 1) 
fi = 2.331 Debye = 2.331 x 10 -18 esu cm 
Ri = (see Section [TOl or. for lio — I 11 , 2n — 2 i 2 , 
or 3 12 — 3 i 3 see Table HD 
gj = 2 J u + 1 

gK = 2 for Ky 0 (using the slightly-asymmetric 
rotor approximation), 

3 

gi = — for K -1 odd 

We can derive the following equation for the molecular 
column density in H 2 CO as derived from a measurement 
of an ortho-H^CO (K_i odd) K-doublet transition as¬ 
suming optically thin emission using Equation 1801 


N tot {H 2 CO) = 


3 h J u ( J u + 1) Q rot 


87 r 3 g, 2 Ri K 2 

ex p (Jl) 


exp 


( kT ex ) 


- 1 


9j9k9i 
f T R dv 


f {J v {T ex ) - MT bg )) 


9.83 X 10 11 J u (Ju + 1 )Qrot 


R,X 2 (2J, l 
ex P (fcfcy 


1) 


f Tndv(km/s) 


exp ( 73 ^) —lL/ {Jv(T e x) Ju(T bg )) 

( 100 ) 


Note that the K-doublet transitions of H 2 CO are rather 
unusual in that, due to an unusual collisional excitation 
effect, these transitions are often measured in absorption 
against the cosmic microwave background radiation. For 
n(H 2 ) < 10 5 ’ 5 cm~ 3 , a collisional selection effect overpop¬ 
ulates the lower ener gy states of the li n — I 11 through 
5-14 — 5i5 transitions ([Evans et al.l 11975t lOarrison et al.1 
I1975H . This overpopulation of the lower energy states 
results in a “cooling” of the J< 5 K-doublets to an exci¬ 
tation temperature that is lower than that of the cosmic 
microwave background. This causes the J < 5 K-doublet 
transitions to appear in absorption. For volume densities 
n(H 2 ) > 10 5 ' 5 cm -3 and Tk — 40 K this collisional pump 
is quenched. For these higher volume densities the J < 
5 K-doublets are driven into emission over a wide range 
of kinetic temperatu res and abundances (see Figure 1 in 
iMangum et al.lf2008f) . 


cm 


16. SUMMARY 

The most general form of the molecular spectral line 
column density is given by Equation 1321 


N tnt = 


3 h Qrot 


87T 3 | im v 


9* 
exp 


■ exp 


E u 

kT P: 


hv 

kT PP 


- 1 


J r v dv. ( 101 ) 


For absorption line measurements one should use the ve¬ 
locity integrated opacity directly in Equation [32j For 
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emission lines, the velocity integrated opacity is usually 
converted to an integrated intensity. In the optically thin 
limit, the column density becomes (Equation 1801) : 


■\fthin _ 

iy tot ~ 


f 3 h 


\8n 3 S [i 2 R. 


/ Qrot \ 
i) \9 j9k9i) 


exp 


\ kT ex J 


9J 9K9i) exp(^)-l 

[ T R dv 


1 


{Ju{T ex ) - MT bg )) J f 


( 102 ) 


The optically thick column density is related to the op¬ 
tically thin column density by Equation 1861 


Ntot 


_ 7 \jthin 

— ly tot 


T 

1 — exp(—r) 


(103) 


In this tutorial, we have carefully derived the terms in 
Equations [3^1 l80l andlSSl Note, though, that these equa¬ 
tions assume all energy levels are characterized by the 
same excitation temperature (the “CTEX LTE approx¬ 
imation” ; See Section fl2l) ). We have also explored the 
limits of the CTEX approximation fSection fl2l) . Non- 
LTE determinations of the column density require a cou¬ 
pled radiative transfer and statistical equilibrium calcu¬ 
lation and can be performed with public.ally available 
tools. 
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APPENDIX 


LINE PROFILE FUNCTIONS 
For a Gaussian profile the function (j){y) is given by 


<t>(y) = 


l 


\[2/ko 


■ exp 


(u- Uq) 

2a 2 


21 


(Al) 


where v$ and v are the central and offset frequencies of the line profile for a particle with mass M and velocity v , 


2er 2 = 

c 2 


2kTi. 


M 


,2 


(A2) 
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and 


J <j>{v)dv = 1. 


The Gaussian profile has a FWHM given by (in both frequency and velocity): 

-i 2 


A VD = 


2 Vi 


A vd = 2 


- In 2 

In 2 


/ 2fcTfc 

^ M 


/ 2fcTfc 

V M 


(A3) 

(A4) 

(A5) 


and a peak value given by: 


fyiy^peak — 


2 v / hT2 
spkAv D 
2-\/ln 2c 
y/TTV 0 Av D ' 


(A6) 


If one uses peak values instead of integrating over a Gaussian profile to derive column densities, one must make the 
following correction: 


In 2 


(Aj 0 j )Gauss — ^ y tot)peak- 

INTEGRATED FLUXES VERSUS BRIGHTNESS TEMPERATURES 


(A7) 


All calculations in this document assume the use of integrated brightness temperatures (/ T B Av). If one uses 
integrated fluxes (/ S v Av), the total molecular column density assuming optically-thin emission ('Equation l82l) is 
modified by using the relationship between flux density and brightness temperature: 


and becomes 


N, nf = 


D 

2R, 


2kT B v 2 

OlJ o ' M 


3c 2 / Q r 


(Bl) 


167T 3 Sfj. 2 v 3 \gjgK9nud e 


exp ( [ S„Av 


kT 


3c 2 


Qrot 


E n 


6XP ' fcT * ' S, ' AV ‘ 


167r 3 f2 s 5/r 2 ^ 3 \gjgKgnud e 

INTEGRATED INTENSITY UNCERTAINTY 


(B2) 


For cases where you do not have a calculation from a fit to the integrated intensity of a spectral line, one can use 
the following estimate given a measurement of the baseline RMS and line profile properties. 


N 


Tdv = Av c ^2 T n 


(Cl) 


n =1 


=/. 


where Av c is the spectral velocity channel width, T n is a spectral channel value, and the line spans N channels. The 
statistical uncertainty of the integrated line intensity is given by: 


N 


!= E' 


n=l 

N 


dl 


dT 

= E a Ty(Aiy ; ) 


2 

a Av 


dl 


d(Av c ) 


n— 1 


= Ncr? (Av c f 


(C2) 


where we have used the fact that we know the velocity channel width (<r Av c = 0) and assuming that all of the channel 
noise values are equal. We can then write Equation IC2I as follows: 


a i = VN cjtAvc 

= a T \J Avu ne Av c 


(C3) 
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where we have used the fact that the spectral line width Avn ne = NAv c to get the last expression for ay. 

HYPERFINE RATIOS AND OPTICAL DEPTH 

For molecules with measurable hyperfine splitting one can derive an estimate of the optical depth in a given transition. 
Starting with the form of the radiative transfer equation which involves the observable Source Radiation Temperature 
Tr ('Equation I27D . we can write the ratio of two hyperfine transitions as follows: 


Ffi(ra) _ f m [J v (T ex<m ) - J u (T bg )] [1 - exp(-T„, m )] 
Tr(s) f s [J v {T eX}S ) - J v (T bg )\ [1 - exp(—r„ iS )] 


(Dl) 


where “m” and “s” refer to the main and satellite hyperfine components of a given rotational transit iorF^I. If we further 
assume that the hyperfine levels are in LTE (Jv( T ex , m ) = J v (T ex ,s)B and that the two hyperfine components emit 
over the same volume of gas (f rn = f s ), Equation IDll becomes: 


T R {m) _ 1 - exp(-7y >m ) 
T r (s) 1 - exp(-T„ jS ) ' 


(D2) 


In equilibrium the optical depths of the hyperfine transitions within a given rotational transition should be proportional 
to their relative intensities: 

lAA = (D3) 

rii,m 

Inserting Equation ID3I into Equation ID2I results in a simple equation which allows for the calculation of the optical 
depth of a hyperfine transition (and, therefore, the rotational transition under study) using just the measured intensities 
of two of its component hyperfine transitions: 


TR{m) 

Tr{s) 


1 - exp(—T„ >m ) 

Ri,, 


1 exp( R i ’ m 'Ci/ ) m) 

Taking the strongest and second-strongest transitions of the NH 3 (1,1) transition, from Table fill we find that: 


Tr(iti) 

Tr(s) 


(NH 3 ( 1,1)) = 


1 - exp(-T„, m ) 


1 — exp(— 


18 


O' 


(D4) 


(D5) 


EXCITATION AND KINETIC TEMPERATURE 

This section is drawn from Appendix A of lMangum et all (119921) . If the metastable (J = K; A t j < 10 - 7 s -1 ) states 
in NH 3 are coupled only through collisions and the populations in the upper non-metastable (J ^ K\ A tJ ~ 0.01s -1 ) 
states in each K-ladder can be neglected, the populations in the metastable states are related through the Boltzmann 
equation. In molecular clouds, though, A K = 1 collisions across K-ladders will deplete metastable states in favor of 
their next lower J metastable states. Therefore, for example, collisional de-excitation of the (2,2) transition will result 
in an increase in the population of the ( 2 , 1 ) state, followed by quick radiative relaxation of the ( 2 , 1 ) state into the 
(1,1) state. This implies that an excitation temperature, T ex (J', K'; J,K) relating the populations in the (J',K') and 
(J,K) states, n(J',K') and n(J,K), may be derived. From the Boltzmann equation 


n(J',K‘) _ g{J', K') 
n(J, K) g(J,K) 


■ exp 


A E(J',K'-J,K) 


T ex ( J', K '\,/, K) J ’ 


(El) 


and the ratio of level (J / ,K / ) and (J,K) column densities (assuming hv -C kT ex (J', K '; J,K )) for the (J,K) and (J'jK 7 ) 
transitions 

N(J\K') J'(J' + 1)K 2 t(J’,K')Av(J',K') 

N(J,K) J(J + 1)(K') 2 t(J, K)Av(J, K) ’ 1 ’ 

and the fact that in a homogeneous molecular cloud 

n(J', K') N(J',K') 


we find that 




■ exp 


g(J,K) 

12 Equivalently, one can use molecular isotopologues with the 
additional assumptions that one knows the isotopic ratio and there 
is no fractionation between that atomic ratio and the molecular 
ratio. 

13 But note that anomalous hyp erfin e excitation has been ob¬ 
served in some molecules (see Section HoTl . notably the HCN J=l—0 


N(J,K) ’ 

J'(J' + 1)K 2 t(J', K')Av(J',K') 
T ex {J ', K'\ J, K) J ~ J(J +1)(K'Yt(J 1 K)Av{J 1 K) ' 


n(J, K ) 
A E(J',K'-J,K) 


(E3) 


(E4) 


transition dLoughnane et ahl 120121) . This would invalidate the 
equal-excitation temperature assumption. 
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Using 




F,F' F,F' 

yy rf,f'/ yy Rm 


T(J,K,m), 


(E5) 


where R is the relative intensity for a quadrupole (F,Fi) or main (to) hyperfine component and 


IjK = 


F,F' f,f' 

yy rf , f > / yy -Rm 


1.0 

-TTI- 

12 ~ 12 

= 2.000 for the (1,1) transition (see Tables fill or fl9l). 
1.0 

_ 56 | 25 j 3 ’ 

135 108 20 

= 1.256 for the (2,K) transitions (see Tables H2l or l20l). 

_ 54 
_ 43’ 

1.0 

_ 45 | 121 , 10 1 

112 ' 432 ' 189 

_ 216 
_ 193 ’ 

= 1.119 for the (3,K) transitions (see Table fl3l). 

1.0 

— 968 | 361 j 35 ’ 

2475 1200 144 

_ 2200 
~ 2057’ 

= 1.070 for the (4,K) transitions (see Table [HI). 


for the (1,1) and (2,2) transitions, we can relate the total optical depth r(J,K) to the op tical depth in the main 
hyperfine component r(J,K,m), noting that for NH 3 g(J,K) = 2 J u + 1, and solving Equation IE4I for T ex (J ', AT'; J, K ) 
we find that 


Using 


T ex (J',K'- J,K) 


-A E(J',K';J,K) 

w f ln [ ( 2 J+ 1)J'(J' + 1)K 2 Ij>k't(J', K', m)Av(J r , K') 
X \ [ (2 J' + 1) J (J + 1)(K') 2 I jk t(J, K, to)A v(J, K) 



TB(J,K,m) 1 — exp [— r( J, K, m)} 
Tb{J',K' , to) 1 — exp [—r( J', A"', to)] ’ 


(E 6 ) 


(E7) 


which assumes equal excitation temperatures and beam filling factors in the (J,K) and (J^K 7 ) transitions. Solving 
Equation IE7I for t(J', K’,m) yields 


r( J 7 , AT', to) = — In 



T B (J', AT', to) 
T b (J, AT, to) 


{l-exp[-T(J,A',TO)]} 


Substituting Equation IE 8 I into Equation IE 6 I 


(E 8 ) 


T ex (J', K'\ J, K) = -A E{J', K'\ J, K) 

[ [ (2J + 1)K 2 Ij, k ,Av{J',K') f 

X [ n [{2J' + 1)J{J + l){K') 2 Ij K T{J,K,m)Av{J,K) V 

For (J',K') = (2,2) and (J,K) = (1,1), Equation IE9I becomes: 


T B (J',K',m) 
Tb(J, AT, to) 


{1-exp [—r( J, A', to)]} 



(E9) 


1^(2,2; 1,1) = -41.5 


In 


0.283Ac(2,2) 
r(l,l,TO)Au(l,l) ln 



Tb{ 2 , 2 , to) 

T b { 1,1 ,m) 


{l-exp[-r(l,l, 




(E10) 
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Table 11 

J-Ijv Hyperfine Intensities for NH3(1,1) 


F[ —> Fi a 

J' ->• J a 

a 

b 

c 

(2U + l)(2Fi + l) 
(21jv+l) 

6j 

Ri (Fi, J) 

(0,1) 

(2d) 

(2,2) 

(1,1) 

(1,1) 

(1,1) 

1 

1 

1 

1 

1 

2 

1 

2 

1 

1 

5 

25 

3 

1 

1 3 

6 i 

1 

§ 

¥ 

12 

(1,1) 

(1,1) 

1 

1 

i 

3 

1 

6 

1 

¥ 

¥ 

9 

(1,2) 

(1,0) 

(1,1) 

(1,1) 

1 

1 

1 

1 

2 

1 

5 

1 

1 

6 1 

3 

a In = I- 












Table 12 




J-Ijv 

Hyperfine Intensities for NH3(2,K) 


F' -► FU 

J' -> J a 

a 

b 

c 

(2F' + l)(2Fi+l) 
(2/w + l) 

6j 

RiUiV) 

(1,2) 

(2,2) 

1 

2 

2 

5 

1 

f° 

15 

1 

¥ 

135 

(3,2) 

(2,2) 

1 

2 

3 

35 

3 

(3,3) 

(2,2) 

1 

3 

2 

49 

3 

2 y/2 

3\/35 

56 

135 

(2,2) 

(1,1) 

(2,2) 

(2,2) 

1 

1 

2 

1 

2 

2 

25 

3 

3 

1 

6 1 

2v/5 

25 

i e 8 

20 

(2,3) 

(2,2) 

1 

2 

3 

35 

3 

1 

7 

1J5 

20 

(2,1) 

(2,2) 

1 

2 

2 

5 

1 

10 

a In = 1. 












Table 13 




J-Ijv Hyperfine Intensities for NH3(3,K) 


-> Fi a 

j’ ->• j a 

a 

b 

c 

(2F 1 '+l)(2F 1 + l) 

(21JV + 1) 

6j 


(2,3) 

(3,3) 

1 

3 

3 

35 

3 

1 

? 1 

28 

5 

189 

(4,3) 

(3,3) 

1 

3 

4 

21 

3 

112 

(4,4) 

(3,3) 

1 

4 

3 

27 

7s 

4^21 

84 

45 

112 

(3,3) 

(3,3) 

1 

3 

3 

49 

3 

121 

432 

(2,2) 

(3,3) 

1 

2 

3 

25 

3 

272 

3735 

40 

189 

(3,4) 

(3,2) 

(3,3) 

(3,3) 

1 

1 

3 

3 

4 

3 

21 

35 

3 

21 

3 

1 U 

189 


a In = 1. 


To derive the gas kinetic temperature from T ex ( 2, 2; 1,1), one uses statistical equilibrium (noting that only collisional 
processes arc allowed between the different K-ladders), detailed balance, and the Boltzmann equation to calculate T k 
from T ex (J', K'\ J, K ). Assuming that the populations in the (1,1) and (2,2) transitions are much greater than that in 
the higher lying levels of para-NH 3 and that the population of the non-metastablc (2,1) level is negligible in comparison 
to that in the ( 1 , 1 ) level, we can use this “three-level model” of NH 3 to analytically derive an expression relating 
T ex { 2 , 2 ; 1 , 1 ) and T K 


1 , <7(2,2; 2,1) 

C(2, 2; 1,1) 


f g( 1,1) AE(2, 2; 1,1) 

15(2, 2) ° XP [ T ex (2, 2; 1,1) 


ff( 2,2) 
<?( 1 , 1 ) 


exp 


A-B(2, 2; 1 ,1) 


T k 


(Ell) 


wher e C(J ', K'\ J, K) is the collisional excitation rate at temperature Tk between levels (J',K') and (J,K). Equation 
IE11I can be re-written as 


T ex ( 2 , 2 ; 1 , 1 ) 




\ | g(2,2;2,l) 
C(2, 2; 1,1) 


T k = 0 . 


(E12) 


Solutions of Equation IE12I give Tk for a measured T ex { 2, 2; 1,1). 


NH 3 FREQUENCY AND RELATIVE INTENSITY CALCULATION TABLES 
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Table 14 

J-Ijv Hyperfine Intensities for NH3(4,K) 


n 


Fi a J' -4- J a 


u „ (2F^+1)(2F 1 +1) 

b C (2/jy+l) 


6j Ri(F 1: J) 


(3,4) 

(4,4) 

1 

4 

4 

21 

1 

? 6 

fi 

180 

7 

IS* 

(5,4) 

(4,4) 

1 

4 

5 

33 

(4,4) 

(4,4) 

1 

4 

4 

27 

361 

1200 

(5,5) 

(4,4) 

1 

5 

4 

121 

3 

2+2 

5+33 

968 

2475 

(3,3) 

(4,4) 

1 

3 

4 

49 

3 

+5 

4yTT 

35 

144 

(4,5) 

(4,4) 

1 

4 

5 

33 

1 

36 

11 

6 £5 

432 

(4,3) 

(4,4) 

1 

4 

4 

21 


In = 1. 


Table 15 

Fi— 1 h Hyperfine Frequencies and Intensities for NH3(1,1) 


F' -4- F a 

F[ -4 Fi a 

/±. vhf ° (kHz) 

a 

b 

C 

(2F' + 1)(2F+1) 
(2/ft+D 

6j 

Ri(F, Fi) 

(I 

(0,1) 

-1568.487 

1 

2 

1 

2 

l 

2 

1 

V6 

1 

3 

(~ ~) 

V2’ 2 ' 

(0,1) 

-1526.950 

1 

2 

i 

3 

2 

4 

1 

/6 

2 

3 

(3 IN 

(2,1) 

-623.306 

1 

2 

2 

3 

2 

4 

2 

V3 

1 

3 

(5 3 \ 

^ 2 ’ 2 ' 

(2,1) 

-590.338 

1 

2 

2 

5 

2 

12 

1 

2\/5 

3 

5 

V2’ 2 ' 

(2,1) 

-580.921 

1 

2 

3 

2 

2 

8 

1 

2v / 30 

1 

15 

(i i) 

!:> 
(1-1) 

(1,1) 

(1,1) 

-36.536 

-25.538 

1 

? 

1 

2 

1 

2 

1 

1 

3 

2 

4 

1 

? 

6 

2 

? 

? 

25 

(2,2) 

-24.394 

2 

2 

12 

1 

10x/3 

(3 3 \ 

V2’ 2 ’ 

(2,2) 

-14.977 

1 

2 

3 

2 

2 

8 

3 

10/2 

18 

50 

''2’ 2 ' 

(1,1) 

+5.848 

1 

2 

i 

3 

2 

4 

"v! 

15. 

1 

9 

(- -) 

' 2 ’ 2' 

(2,2) 

+10.515 

1 

2 

5 

2 

2 

18 

14 

25 

(3 3 \ 

V2’ 2 ' 

(1,1) 

+16.847 

1 

2 

3 

2 

1 

8 

+5 

6+2 

5 

9 

V2’ 2 ' 

(2,2) 

+19.932 

1 

2 

2 

5 

2 

12 

1 

10\/3 

1 

25 

^ 2 ’ 2' 

(1,2) 

+571.792 

1 

2 

2 

3 

2 

4 

1 

2^3 

2+30 

1 

3 

(3 3 \ 
v 2 ’ 2 / 

(1,2) 

+582.790 

1 

2 

3 

2 

2 

8 

1 

15 

f- -) 

v 2 ’ 2 > 

(1,2) 

+617.700 

1 

2 

2 

5 

2 

12 

1 

2\/5 

3 

5 

V 2 ’ 2 / 

(1,0) 

+ 1534.050 

1 

2 

1 

2 

i 

2 

1 

%/6 

1 

3 

(3 In 
^ 2 ’ 2' 

(1,0) 

+ 1545.049 

1 

2 

i 

3 

2 

4 

1 

n/6 

2 

3 


a F = Ih + Fi and F\ = J + In, where In = 1 and Ih = J. 

b Frequency offset in kHz relative to 23694.495487 kHz (IKukolicblll967l . Table I (cal¬ 
culated)). 
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Table 16 

Fi— 1 h Hyperfine Frequencies and Intensities for NH3(2,2) 


F' -+ F a 

F[ ->• Fi a 

A vhf° (kHz) 

a 

b 

C 

(2F'+1)(2F+1) 

(2/ft+II 

6j 

Ri(CFi) 

(3 3 \ 

V2’ 2' 

(1,2) 

-2099.027 

1 

2 

3 

2 

2 

8 

1 

2\/30 

1 

15 

(- 

V 2 5 2/ 

(1.2) 

-2058.265 

1 

2 

2 

5 

2 

12 

1 

2v/5 

3 

5 

(- ~) 

v 2 ’ 2 t 

(1.2) 

-2053.464 

1 

2 

2 

3 

2 

4 

1 

2v^3 

1 

3 

(1 

^ 2 ’ 2 ' 

(3,2) 

-1297.079 

1 

2 

3 

7 

2 

24 

1 

^42 

+53 

4 

7 

(5 3 \ 

V2’ 2 > 

(3,2) 

-1296.096 

1 

2 

3 

5 

2 

12 

6 

15 

(if) 

(3,2) 

-1255.335 

1 

2 

5 

2 

3 

18 

1 

Sv'TO 

1 

35 

(§,§) 

(I, I) 

(1,1) 

-44.511 

1 

1 

2 

1 

3 

4 

1 

6 

1 

? 

25 

(2,2) 

-41.813 

2 

2 

12 

1 

10v/3 

(- -) 

V2’ 2' 

(3,3) 

-41.444 

1 

2 

3 

7 

2 

24 

1 

14\/6 

1 

49 

(§. 

ri 

k 2 ’ 2 ' 

(2,2) 

(2,2) 

-1.051 

-1.051 

1 

1 

2 

5 

2 

2 

2 

18 

8 

+7 

10+2 

14 

¥ 

25 

(I 

k 2 ’ 2' 

(3,3) 

+0.309 

1 

2 

7 

2 

3 

32 

3+3 

28+2 

27 

49 

(- -) 

v 2 ’ 2 > 

(3,3) 

+0.309 

1 

2 

5 

2 

3 

18 

V 10 

21 

20 

49 

(if) 

(1,1) 

+1.054 

1 

2 

3 

2 

1 

8 

>/5 

6+2 

5 

9 

(if) 

(if) 

(1,1) 

(2,2) 

+1.054 

+39.710 

1 

1 

2 

1 

2 

2 

1 

5 

2 

2 

12 

1 

_L 

10+3 

2 

? 

25 

(if) 

(3,3) 

+42.045 

1 

2 

3 

7 

2 

24 

1 

14v/6 

1 

49 

(if) 

(■2’ 2 ) 

(1,1) 

(2,3) 

+46.614 
+ 1254.584 

1 

! 

2 

1 

5 

2 

3 

2 

3 

4 

18 

1 

1 5 

3+70 

1 

? 

35 

(3 5x 

v 2 ’ 2 / 

(2,3) 

+ 1295.345 

1 

2 

3 

5 

2 

12 

1 

v/30 

6 

15 

(5 7x 
v 2 ’ 2 l 

(2,3) 

+ 1296.328 

1 

2 

3 

7 

2 

24 

1 

\/42 

2+3 

4 

7 

(3 1\ 
v 2 5 2 / 

(2,1) 

+2053.464 

1 

2 

2 

3 

2 

4 

1 

3 

(5 3 \ 

V 2 ’ 2 / 

(2,1) 

+2058.265 

1 

2 

2 

5 

2 

12 

1 

2>/5 

3 

5 

(3 3 \ 

^ 2 ’ 2 ' 

(2,1) 

+2099.027 

1 

2 

3 

2 

2 

8 

1 

2v / 30 

1 

15 


a F = Ih + F\ and Fi = J + In, where In = 1 and Ih = h- 

b Frequency offset in kHz relative to 23722633.335 kHz iKukolicbl 119671 . Table II 


(calculated)) 
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Table 17 

Fi— 1h Hyperfine Frequencies and Intensities for NFI 3 (3,3) 


F' —> F a 

F' ->■ Fi a 

AisHF h (kHz) 

a 

b 

c 

(2F' + 1)(2F+1) 
(2-Tff + l) 

6j 

Ri(F,Fi) 

(1 9) 

V2 ’ 2' 

(2,3) 

-2324.577 

3 

2 

3 

9 

2 

20 

1 

2v / 14 

5 

14 

(- -) 
ri 

V2 ’ 2' 

(2,3) 

-2312.558 

3 

3 

7 

12 

1 

12 

49 

(2,3) 

-2304.415 

2 

3 

2 

2 

1 

2v / 5 

1 

10 

(3 5n 

V 2 ’ 2 ' 

(2,3) 

-2301.989 

3 

2 

3 

5 

2 

6 

+2 

5+3 

4 

25 

(ID 

(2,3) 


3 

2 

7 

2 

3 

16 

1 

14\/2 

2 

49 

(f.f) 

if) 

(14) 

(2,3) 

(2,3) 

(2,3) 


3 

2 

5 

2 

3 

3 

3 

9 

4 

12 

8 

J05 

¥ 

14+30 

64 

12^5 

¥ 

490 

(f.f) 

(2,3) 


3 

2 

3 

2 

3 

6 

1 

10+51 

1 

350 

(- -) 
v2 ’ 2' 

(4,3) 

-1690.763 

3 

2 

4 

7 

2 

12 

5 

28\/2 

150 

784 

(9 7n 

V 2 ’ 2 ' 

(4,3) 

-1689.154 

3 

2 

4 

9 

2 

20 

vTl 

12+6 

110 

432 

(- -) 

V2 ’ 2' 

(4,3) 

-1682.925 

3 

2 

4 

5 

2 

6 

1 

\/42 

1 

7 

(¥4) 

(4,3) 

-1679.029 

3 

2 

4 

11 

2 

30 

1 

3v/T0 

1 

3 

(f.f) 

(4,3) 


3 

2 

9 

2 

4 

25 

1 

6\/30 

5 

216 

(1 1) 

V2’ 2' 

(4,3) 


3 

2 

7 

2 

4 

16 

+5 

21+6 

40 

1323 

(5 5 \ 

V2 » 2' 

(4,3) 


3 

2 

5 

2 

4 

9 

1 

14\/2 

9 

392 

(5 7 n 

(4,3) 


3 

2 

3 

7 

2 

12 

1 

84^2 

1 

1176 

(1 9 N 

V2 ’ 2^ 

(4,3) 


3 

3 

9 

20 

1 

1 


2 

2 

12+210 

1512 

(- -) 

V2 ’ 2' 

(2,2) 


3 

2 

2 

7 

2 

12 

1 

v / 2l0 

6 

105 

(5 3 n 
n 2 ’ 2 ' 

(2,2) 

-80.030 

3 

2 

2 

5 

2 

6 

+7 

10+6 

7 

100 

(3 In 

V2 ’ 2' 

(2,2) 


3 

2 

2 

3 

2 

2 

1 

2VI0 

1 

20 

(f.f) 

(3,3) 

-64.182 

3 

2 

3 

5 

2 

6 

1 

v/210 

3 

105 

(f.f) 

(3,3) 


3 

2 

3 

7 

2 

12 

+5 

28+2 

15 

392 

(9 7 n 

(3,3) 

-61.951 

3 

2 

3 

9 

2 

20 

1 

4\/42 

5 

168 

fii 9 N 

v 2 ’ 2' 

(4,4) 

-50.048 

3 

2 

4 

11 

2 

30 

1 

5\/66 

1 

55 

(9 7 n 

V 2 ’ 2 ' 

(4,4) 


3 

2 

4 

9 

2 

20 

+77 

180-^2 

_ 4+42 

77 

3240 

(I 5n 

V 2 ’ 2 ' 

(4,4) 


3 

2 

4 

7 

2 

12 

1 

56 

(f.f) 

(3,3) 

-1.334 

3 

2 

3 

2 

3 

4 

1 

\/35 

4 

35 

(f.f) 

(3,3) 


3 

2 

9 

2 

3 

25 

+TT 

2+210 

275 

840 

(1 1) 

V2’ 2' 

(3,3) 


3 

2 

7 

2 

3 

16 

n/2 

42+10 

32 

147 

(5 5n 

V2 ’ 2' 

(3,3) 


3 

2 

5 

2 

3 

9 

2601 

17640 

(- -) 

' 2 ’ 2 ' 

(4,4) 

+0.445 

3 

2 

5 

2 

4 

9 

_ 5 

6\/42 

25 

168 

(f.f) 

(4,4) 


3 

2 

7 

2 

4 

16 

4+2 

9+35 

512 

2835 

(f.f) 

(4,4) 


3 

2 

9 

2 

4 

25 

41 

90+22 

1681 

7128 

(ii I!) 
v 2 ’ 2 ^ 

(4,4) 


3 

2 

11 

2 

4 

36 

+13 

3+165 

52 

165 

(I I) 
v 2 ’ 2 ' 

(2,2) 

+1.067 

3 

2 

1 

2 

2 

1 

1 

2^5 

1 

20 

(3 3) 

V 2 ’ 2 ' 

(2,2) 


3 

2 

3 

2 

2 

4 

1 

5\/2 

2 

25 

(5 5n 

V2 ’ 2' 

(2,2) 


3 

2 

5 

2 

2 

9 

11 

3y/7 

1089 

6300 

(1 

'2 ’ 2' 

(2,2) 


3 

2 

7 

2 

2 

16 

+3 

2+35 

12 

35 

(£ 11) 
V2’ 2 ' 

(4,4) 

+50.937 

3 

2 

4 

11 

2 

30 

1 

1 

55 

(I 9) 

' 2 ’ 2 ' 

(4,4) 


3 

2 

4 

9 

2 

20 

+77 

180^/2 

4+42 

77 

3240 

(5 7 n 

V2 ’ 2' 

(4,4) 


3 

2 

4 

7 

2 

12 

1 

56 

(1 9) 
v2 ’ 2' 

(3,3) 

+59.283 

3 

2 

3 

9 

2 

20 

1 

4v/42 

5 

168 

(5 7 n 
v 2 ’ 2 ' 

(3,3) 


3 

2 

3 

7 

2 

12 

V5 

28n/2 

15 

392 

(3 5n 
n 2 ’ 2 ' 

(3,3) 


3 

2 

3 

5 

2 

6 

1 

\/2l0 

3 

105 

(I 3) 

(2 ’ 2' 

(2,2) 


3 

2 

2 

3 

2 

2 

_ 1 

2 v / 10 

1 

20 

(3 5n 
v 2 ’ 2 ' 

(2,2) 

+82.164 

3 

2 

2 

5 

2 

6 

+7 

10+6 

7 

100 

(5 7 n 
^2 ’ 2' 

(2,2) 


3 

2 

2 

7 

2 

12 

1 

\/2l0 

6 

105 



Continued 

on 

next 

page 
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Table 17 

Fi-I# Hyperfine Frequencies and Intensities for NH3(3,3) (continued) 


F' -+ F a 

-+ Fi a 

Ai y HF b (kHz) 

a 

b 

C 

(2F' + 1)(2F+1) 
(21rr+l) 

6j 

Ri(CFi) 

(if) 

(3,4) 

+ 1682.036 

3 

2 

4 

5 

2 

6 

1 

■s/42 

1 

7 

(if) 

(3,4) 

+1688.264 

3 

2 

4 

9 

2 

20 

V+T 

12+6 

55 

216 

(§>£) 

(3,4) 

+1689.873 

3 

2 

4 

7 

2 

12 

_ 5 

2S>y/2 

75 

392 

V 2 ’ 2 > 

(3,4) 

+ 1678.140 

3 

2 

4 

11 

2 

30 

1 

3-\/T0 

1 

3 

(9 I) 
v 2 ’ 2 ' 

(3,4) 


3 

7 


20 

1 

1 


2 

2 


12+210 

1512 

(I 

V2 ’ 2/ 

(3,4) 


3 

2 

5 

2 

4 

12 

1 

84\/2 

1 

1176 

9 \ 

V2 ’ 2/ 

(3,4) 


3 

2 

9 

2 

4 

25 

1 

6v / 30 

5 

216 

(I 1) 

V 2 ’ 2 / 

(3,4) 


3 

2 

7 

2 

4 

16 

+5 

21+6 

40 

1323 

f- -) 

v 2 ’ 2 ' 

(3,4) 


3 

2 

5 

2 

4 

9 

1 

14v/2 

9 

392 

(5 3) 

V2 ’ 2' 

(3,2) 

+2301.723 

3 

2 

3 

5 

2 

6 

\/2 

5n/3 

4 

25 

V2 ’ 2' 

(3,2) 

+2304.148 

3 

2 

3 

3 

2 

2 

i 

2^5 

2+14 

1 

10 

(I 5x 

r? 

V2 ’ 2’ 

(3,2) 

(3,2) 

+2312.291 

+2324.310 

3 

2 

3 

3 

7 

2 

12 

20 

12 

f 

14 

(k 1) 

V2 ’ 2' 

(3,2) 


3 

2 

2 

7 

2 

12 

1 

14+30 

1 

490 

(f.f) 

(3,2) 


3 

2 

2 

5 

2 

6 

1 

10^21 

j.05 

14V2 

1 

350 

(f,f) 

[?’?) 

\2 ’ 2/ 

(3,2) 

(3,2) 

(3,2) 


3 

| 

2 

3 

? 

2 

3 

3 

3 

4 

9 

16 

i 

!^5 

49 


a F = Ih + Fi and Fi = J + In, where In = 1 and Ih = f- 

b Frequency offset in kHz relative to 23870129.183 kHz (IKukolichll 19671 . Table III (cal¬ 


culated)). 
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Table 18 

Fi— 1h Hyperfine Frequencies and Intensities for NH3(4,4) 


F' _». F a 

F( Fi a 

C\v HF h (kHz) 

a 

b 

C 

(2F' + l)(2f+l) 
(2+r+l) 

6j 

Ri(F,Fi) 

(1 9) 

V2 » 2' 

(3,4) 

-2485.883 

1 

2 

4 

9 

2 

40 

_1 

6\/2 

5 

9 

(- ~) 

V2 ’ 2/ 

(3,4) 

-2481.291 

1 

2 

4 

7 

2 

24 

1 

2y/l4 

3 

7 

C- 

\ 2 ’ 2 ' 

(3,4) 


1 

2 

7 

2 

4 

32 

1 

12^+4 

3+TO 

1 

63 

(9 7x 

V2 ’ 2' 

(5,4) 

-1924.144 

1 

2 

5 

9 

2 

40 

4 

9 

V 2 ’ 2' 

(5,4) 

-1922.222 

1 

2 

5 

11 

2 

60 

1 

\/Tio 

6 

11 

(£ 9x 
\2’ 2' 

(5,4) 


1 

2 

9 

2 

5 

50 

1 

15+22 

1 

99 

f9 9x 

(4,4) 


1 

2 

9 

2 

4 

50 

VT2 

45 

44 

81 

(1 I) 

V2 ’ 2' 

(4,4) 

-1.448 

1 

2 

7 

2 

4 

32 

+35 

36+2 

35 

81 

(9 

V2 ’ 2' 

(4,4) 


1 


9 

40 

1 

1 


2 


2 

lsyTo 

18+10 

81 

(I £) 

V2 ’ 2' 

(4,4) 


1 

2 

4 

9 

2 

40 

1 

81 

(ii I!) 

(5,5) 


1 

2 

11 

2 

5 

72 

+65 

66+2 

65 

121 

(£ £) 
(A’ 

v 2 ’ 2' 

(5,5) 

(5,5) 

+0.526 

1 

? 

9 

2 

5 

n 

50 

60 

3+3 

¥ 

54 

n 1 


2 


2 

22-^15 

22v/l5 

3y3 

28\/2 

121 

(9 11) 
v 2 ’ 2 ' 

(5,5) 


1 

2 

5 

11 

2 

60 

1 

121 

(I-D 

(3,3) 


1 

2 

7 

2 

3 

32 

27 

49 

( f-f } 

(I-1) 

(3,3) 

(3,3) 


1 

? 

2 

5 

2 

3 

3 

7 

2 

18 

24 

+To 

¥ 

14+6 

20 

¥ 

49 

V2 ’ 2' 

(3,3) 


1 

2 

3 

7 

2 

24 

1 

14^6 

1 

49 

(9 11) 

^ 2 ’ 2 ' 

(4,5) 

+ 1931.301 

1 

2 

5 

11 

2 

60 

1 

\/TIo 

6 

11 

f ) 

V2 ’ 2/ 

(4,5) 

+ 1923.222 

1 

2 

5 

9 

2 

40 

l 

3\/T0 

4 

9 

(9 9n 

V2 ’ 2/ 

(4,5) 


1 

2 

9 

2 

5 

50 

1 

15-^22 

2+T4 

1 

99 

(Z 5) 

V2 ’ 2' 

(4,3) 

+2480.878 

1 

2 

4 

7 

2 

24 

3 

7 

(9 7) 

V2 ’ 2' 

(4,3) 

+2485.470 

1 

2 

4 

9 

2 

40 

1 

6v/2 

5 

9 

(1 

V2 ’ 2' 

(4,3) 


1 

2 

7 

2 

4 

32 

1 

12+14 

1 

63 


a F = Ih + Fi and Fi = J + In, where In = 1 and Ih = §■ 

b Frequency offset in kHz relative to 24139416.34 kHz (IKukolich fc Wofsvll 19701 . Table II 


(calculated)) 

















How to Calculate Molecular Column Density 


35 


Table 19 

Hyperfine Intensities for NH3(1,1) 


F' -> F F[ —>• F\ J' -> J Ri (Fi, J)Ri (F, Fi ) a ’ b 


(if) (o.i) (i,i) 

(1,1) (o,i) (i,i) 

(4,4) (2,1) (1,1) 

(f,f) (2,1) (1,1) 

(1,1) (2,1) (1,1) 

(1,1) (W) (CD 

(1,1) (M) (CD 

(1,1) (2,2) (1,1) 

(|,f) (2,2) (1,1) 

(1,1) (CD (CD 

( 4 , 4 ) ( 2 , 2 ) ( 1 , 1 ) 

(4,4) (i,i) (i,i) 

(4,4) (2,2) (1,1) 

(4,4) ( 1 , 2 ) (i,i) 

(4,4) ( 1 , 2 ) (i,i) 

(4,4) ( 1 , 2 ) (i,i) 

(4,4) (i,o) (i,i) 

(§, |) ( 1 , 0 ) ( 1 , 1 ) 


¥ 

X 

1J8 

¥_ 

108 

lj 8 

T 98 

1 
1^8 
JF 

198 

JL 

198 

195 

2 

¥ 


a Compare with IKukolichl (119671 ) T able IX aft e r scal - 
ing Ri by (21 H + 1)(2I N + 1) = 6 (IKukolichl (119671 ) 
lists unnormalized line strengths in their Table IX). 
b Note that the sum of the relative intensities 
=* 1 - 0 - 


Table 20 

Hyperfine Intensities for NH3(2,2) 


F' -> F F[ —>• F], J'->J Ri (Fi, J)Ri (F, Fi ) a ' b 


(f, f) (1,2) (2,2) 

(4,4) (1,2) (2,2) 

(4,4) (1,2) (2,2) 

(4,4) (3,2) (2,2) 

(4,4) (3,2) (2,2) 

(4,4) (3,2) (2,2) 

(4,4) (1,1) (2,2) 

(4,4) (2,2) (2,2) 

(4,4) (3,3) (2,2) 

(4,4) (2,2) (2,2) 

(4,4) (2,2) (2,2) 

(4,4) (3,3) (2,2) 

(4,4) (3,3) (2,2) 

(4,4) (1,1) (2,2) 

(4,4) (1,1) (2,2) 

(4,4) (2,2) (2,2) 

(4,4) (3,3) (2,2) 

(4,4) (1,1) (2,2) 

(4,4) (2,3) (2,2) 

(4,4) (2,3) (2,2) 

(4,4) (2,3) (2,2) 

(4,4) (2,1) (2,2) 

(4,4) (2,1) (2,2) 

(f,|) (2,1) (2,2) 


300 

T 90 

W 

6J5 

6|5 


995 


199 

¥ 


995 


195 


190 

300 


a Compare with IKukolichl (119671 ) T able IX aft e r scal - 
ing Ri by (2I ff + l)(2/jy + 1) = 6 (IKukolichl (119671 ) 
lists unnormalized line strengths in their Table IX). 
b Note that the sum of the relative intensities 
J2iRi = CO. 
























